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Dr. Aaron Temkin and Dr. Richard J. Drachman have made significant contributions in the 
field of atomic physics during nearly the last fifty years. It is not possible to enumerate all of their 
contributions but I will describe a few, first those of Temkin and then of Drachman. Some of the 
topics of their research are discussed in the articles that follow this one. 

Method of Polarized Orbitals 

Their work, as is the work of most atomic theorists, is concerned with solving the Schrodinger 
equation accurately for wave function in cases where there is no exact analytical solution. In 
particular, Temkin is associated with electron scattering from atoms (and ions). When he started 
there already were a number of methods to study the scattering of electrons from atoms. The 
simplest approximation for a hydrogenic target is the static approximation 


tf(ri,r 2 ) = w(ri)$target(r 2 ), (1) 

where u( ri) is the scattering function. This gives rise, in a partial wave expansion, to differential 
equations which can be solved easily. An important improvement on this approximation is the 
exchange approximation of Morse and Allis [1] in which the target is assumed to remain in its 
original state in the presence of the incident electron but the ansatz for the wave function is 
explicitly (anti) symmetrized 

^(ri, r 2 ) = u(ri)$target(r 2 ) ± u(r 2 )$target(ri), (2) 

where the upper sign refers to singlet scattering and the lower sign to the triplet scattering. This 
gives rise to an integro-differential equation which will be given later. 

But in reality, the target does not stay in its original state and various forms of polarization 
potentials have been added to the equation for u(r ) to take into account this distortion. The 
method of polarized orbitals proposed by Temkin [2] was the first method to include the effect of 
polarization in the ansatz for the total wave function. It has been used for scattering of slowly 
moving electrons from various targets so that the target is distorted adiabatically. This means the 
target is assumed to follow the instantaneous motion of the scattered electron. The basic problem 
is how to take into account this distortion mathematically. I will discuss it for a simple target like 
the hydrogen atom in its ground state (Is), following Temkin [3]. The Hamiltonian H (in Ry units) 
for two electrons in the field of the nucleus of charge Z is given by 


Jf = - V? - 


2Z_ 

ri 


2 Z 2 

I > 

^2 r 12 


which for an incoming electron 1 can be written as 


( 3 ) 


H = H 0 + V(t u r 2 ), 


( 4 ) 


1 
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where Hq also contains V 2 


TT 2 2 2 Z 2(Z — 1) 

= L . 

r 2 fi 


2 2 

V(ri,r 2 ) = 1 , 

n r 12 


(5) 

( 6 ) 


and ri 2 = |ri — r 2 1 . The first order perturbed wave function of the electron 2 due to the potential 
V can be written as 


^ P ol(ri,r 2 ) = $ioo(r2) - 


/ 

£ 


< n/0|V(ri, r 2 |100 > . . 

2 $n/o(r 2 ), 

€ n ~ ^1 



where integration over continuum states is implied and e n = —Z 2 n 2 is the energy of the nth 
hydrogenic state <f> n i m which is given by 


<Wr 2 ) = (8) 

r 2 

The perturbing potential V for ri > r 2 is given by 

V (i'll r 2 ) = 2 ^2 -j^-Pi(cos0 2 ), (9) 

1= 1 r i 

where we have fixed the z-axis in the direction of iq and have used the well-known expansion of 
l/ri 2 - Noting that the l = 0 term drops out of V(r!,r 2 ), we can write the matrix in Eq. (7) as 


< nlO\V (r x , r^llOO > l/o(^2) = 


The second integral vanishes in the limit r\ — > 00 and the above matrix for large r\ can be 
approximated by 

-E t F / ( cos<? 2 ), ( 11 ) 

(47t)2 

which can be written as 


<n/0|V(n,r 2 )|100>y/ O (fi 2 ) = 


-HT / u n/( r 2 ) r 2 u lo( r 2 )^ r 2 

Jo 


[-nr [ u n i{r' 2 )r'juio{r' 2 )dr' 2 
ry Jo 

f°° /I 11 2 

+ 4 / w n/(r 2 )-^ttio(r 2 )dr 2 ]- — “i" ( cos0 2 ) • (10) 

M r 2 + (47r)5 


< n/0|F( ri ,r 2 )|100 > Y i0 (ft 2 ) = 
Now Eq. (7) can be written as 


L r l 


< n/|r 2 |10 > 


J (4?r) 


-P/(cos0 2 ). 


(12) 
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\Tr /_ rl - * V' 2 P '( cos02 ) V' < nl\r'l\lO > u nl {r 2 ) 

^poi( r i,r 2 ) = $(n, r 2 ) - $100 2 ^ J+T— - L — _ 


To sum this series, use the fact that the bound states satisfy 


V 72 2Z ( \ «n/(r 2 ) D/ 

-V 2 tt n /(r 2 ) = Cn P/(cOS0 2 ), 

r 2 J r 2 


and therefore 


— V2 - ^ - €1 u n i{ r 2 ) = (e„ - t l ) Unl ^ Pi{cosS 2 ), 


which implies (letting A4> = $100 - $(ri,r 2 )) 


= El+i E < n/|r 2 '|10> 


i=l ' 1 n=l+l 

00 2 00 


Unljr 2) F( cos 6> 2 ) 
r 2 (4 tt)^ 


= E~ 7 +T E KK r 2 ) >< nl\r'l \10 > 


1=1 ' 1 n=l + 1 


r 2 (47r) 


Use now completeness 


K/(r 2 ) >< «„;(r 2 )| = £(r 2 - r 2 ), 


for any /, so that A4> satisfies the differential equation 


00 r l 

V—v I O 


-Vi Cl A$ = 2^- ? | r $ 1 oo(r 2 )F ! (cos0 2 ), 

r 2 J 1=1 r l 


where we have used Eq. (8) for 4>ioo- Expanding 


a i _ V - ' ( r 2)/ r 2 F (CPS^ 2 ) 

~ h r i +1 (4 Z7T)I ' 


reduces Eq. (18) to the set of uncoupled equations for each l 


cl 2 2Z 1(1+ 1) 1 . r- 1 , 

— + 2 e i u la -n(r 2 ) = 2VZr 2 u ls (r 2 ), 

ai'2 1 2 r 2 


where ui s (r) = 2^/Z^re Zr . Eq. (20) can be solved analytically 


/ r ^+ 2 

“,»'M= 2 '- Zr2 


3 


Symposium on Atomic & Molecular Physics 


The method of polarized orbitals uses only the dipole part (/ = 1) of A$ [4] so that the polarized 
part of the target is 




pol 


$100 - 


€ ( r ii r 2 ) Uis^ P {r 2 ) Pi(cos6>i 2 ) 


1 1 

12 


n r 2 (AZn)i 

where 6\ 2 is the angle between ri and r 2 and Temkin has introduced the step function 


(22) 


e(ri,r 2 ) = 1, r 1 >r 2 

= 0, ri < r 2 . (23) 

which ensures that polarization takes place only when the scattered electron is outside the orbital 
electron. The second term in Eq. (22) is called $^ ol (ri,r 2 ) and gives the distortion of the orbit of 
hydrogen atom in its ground state (Is). This then solves the problem of the inclusion of polarization 
in the wave function, for substituting Eq. (22) into Eq. (13) one gets the explicit form of the 
polarized orbital wave function 

$(ri, r 2 ) = u/(ri) [$ioo(r 2 ) + $^ ol (n, ^ 2 )] ± «/(r 2 ) [$ioo(ri) + $i ol ( r 2 i ri)] , (24) 

where 


Mr) = (25) 

The scattering equation is obtained from 

jY{ o {Q 1 )*Z(r 2 )(H-E)V(r 1 ,r 2 )d{l 1 dr 2 = 0, (26) 

E=-Z 2 +k 2 , (27) 

where H is the Hamiltonian, E is the total energy, k 2 is the energy of the incident electron and 
$100 = $o. We will derive the scattering equation for any Z, the charge of the nucleus. If the 
variational principle were used it would require that we operate on the left by ($o + ^i o1 )- This 
would give us higher order terms which are not in accord with the first-order perturbation theory. 
(It would also give quadratic terms involving $^ g1 which would be singular if the step function 
e(ri,r 2 ) were retained in its definition.) Therefore, Temkin has multiplied on the left by $ 0 °nly 
in Eq. (26). (In retrospect, the main disadvantage of this ansatz is the loss of any bound on the 
calculated phase shifts from this equation.) The resulting integro-differential equation for tq(r) [5], 
including the p - wave (/=!) correction of Sloan [6] is, letting ri=r, 


d 2 

[ i^ + 


k 2 +v st {r) - — ]^;(r) ± 4 Z 3 e Zr [(k 2 + Z 2 )rSi 0 J dr 2 e Zr2 r 2 ui(r 2 ) 

O pr poo p~Zr2 

~ 2l+l ^ r ~ l Jo dr 2 e~ Z, ' 2 r 2 +l u t (r 2 ) +r l+1 J dr 2 —p—U[(r 2 ))] 


4 



Symposium on Atomic & Molecular Physics 


=F 

+ 


+ 


± 


&L Ul (r) 

(Zry ui[r} 


4 Ze- Zr [\Su{- l -Z{Z 2 + k 2 )r z + (Z 2 - k 2 




dro 


-Zr 2 


-ui(r 2 ) 


l 


2(— + r 2 K — 

^ 2 ^ A (2/-l)(2/ + l) 


/»o 

J -'l 


dr 2 - 


o~ Zr 2 

►J + l 


■u;(r 2 ) 


/+1 ,+i Z 00 , e" Zr » 

(2/ + 1) (2/ + 3) r + J r ^^(r 2 ))] 

^Ze~ 2i Hi[(^Z 2 r 2 + ^Zr - 3)ti/(r) - (^- + r)^] f 


(28) 


where the static potential is given by 


v st (r) = 2^— ^ + 2(Z + -)e~ 2Zr , 
r r 


(29) 


and the direct polarization potential is given by 


f3(Zr)/(Zr ) 4 -> a/r 4 for r -> oc, 


(30) 


where 


/?(*) = \ - | e - 2 "(^ 5 + + 9* 3 + f ^ 2 + f * + f )• ( 31 ) 

Note that the polarized orbital ansatz Eq. (24) provides a natural cut off for the polarization 
potential f3(Zr)/r 4 -» 0 as r -> 0 and gives the correct polarizability as r — >■ oo, where a = 
9/2Z 4 is the dipole polarizability of the target with nuclear charge Z. We get the equation for the 
exchange approximation [1] by putting the right hand side equal to 0 and we get exchange adiabatic 
approximation by retaining the first term on the right hand side which is the direct polarization 
potential. 

Eq. (28) can be solved for the function ui(r) with the phase shifts rji being obtained from the 
values of the function at large distance: 


7 r 

lim r _>ooUi(r) oc sin(kr -\— + rfr) (32) 

The phase shifts for electron-hydrogen scattering for S - wave and P-wave are given in Table I for 
various k in three different methods: exchange approximation, polarized orbital method and Kohn 
variation principle [7,8]. The effect of the polarization is dominated by the direct - a/r 4 potential 
and always has the effect of increasing the phase shifts relative to the exchange approximation. The 
effect of exchange polarization terms is smaller and can be either positive or negative depending on 
the spin. The X S and 3 S results decrease with increasing k while 3 P results increase with increasing 
k. The effect of the exchange polarization terms in the 1 P case leads to three changes in sign of 
phase shift as k increases, indicating that the method does include the essential physics. This is 
the first method to give three changes of sign of the l P phase shifts correctly (cf. Fig. 6-1-5 in [9]). 
We will return to the question of bounds in another section. 
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TABLE I. e-H phase shifts of 1 S', 3 S, 1 P, and 3 P for various k in different approximations. 


Partial Wave 

k 

Exch. 

PO“ 

KVP 6 

Exch. 

PO“ 

KVP 6 



Singlet 



Triplet 



o c 

8.10 

5.9 

5.965 

2.35 

1.9 

1.7686 


0.1 

2.396 

2.583 

2.553 

2.907 

2.945 

2.939 


0.2 

1.870 

2.144 

2.0673 

2.679 

2.732 

2.7171 


0.3 

1.508 

1.750 

1.6964 

2.461 

2.519 

2.4996 

5 

0.4 

1.239 

1.469 

1.4146 

2.257 

2.320 

2.2938 


0.5 

1.031 

1.251 

1.202 

2.070 

2.133 

2.1046 


0.7 d 

0.7415 

0.947 

0.930 

1.748 

1.815 

1.7797 


oo 

O 

0.6513 

0.854 

0.886 

1.617 

1.682 

1.643 


0.1 

-0.0012 

0.0067 

0.007 

0.00220 

0.0109 

0.0114 


0.2 

-0.0084 

0.0171 

0.0147 

0.01666 

0.0486 

0.0450 


0.3 

-0.0240 

0.0210 

0.0170 

0.05192 

0.1151 

0.1063 

P 

0.4 

-0.0461 

0.0163 

0.0100 

0.10497 

0.2005 

0.1872 


0.5 

-0.0703 

0.0064 

-0.0007 

0.16935 

0.2867 

0.2705 


0.6 

-0.0919 

-0.0039 

-0.009 

0.23183 

0.3574 

0.3412 


0.7 

-0.1077 

-0.0100 

-0.013 

0.28329 

0.4063 

0.3927 


0.8 

-0.1154 

-0.0095 

-0.004 e 

0.32044 

0.4351 

0.427 


0.9 

-0.1146 

-0.0015 


0.34436 

0.4487 



1.0 

-0.1058 

0.0135 


0.35791 

0.4520 



a *S-wave phase shifts obtained by Temkin and Lamkin [5] using the method of polarized orbital. 
The P - wave results are from Sloan [6]. 

b Kohn variational results for 5-wave are from Schwartz [7] and P - wave results are from Armstead 

[ 8 ]. 

c The k = 0 entries are scattering lengths [5]. 

^The polarized orbital results are from [9]. 

e Beyond k— 0.8 the phase shift becomes positive and in fact goes through a ( X P) resonance at 
£-0.846 

On the other hand, scattering length a, which is defined as 


limk-m kcot r) = -1/a, (33) 

does have an upper bound, in the variational calculation [7] a exa ct < ^calculated- The scattering 
lengths in various approximations are given in Table I. For the triplet state, Rosenberg et al. [10] 
had obtained a t < 1.91 for the electron-hydrogen scattering. They used an s - wave function having 
a correct asymptotic form at infinity but did not have a slowly vanishing part (^), as subsequently 
pointed out by Temkin [11]. Furthermore, Temkin [11] showed that the expression for the scattering 
length calculated at a finite distance R should be corrected 

z / 1 a 1 a 2 \ , 

o = (34) 


6 


where a is the true scattering length, and a the polarizabilty of the target. With this modification 
Temkin obtained an improved value 
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a t = 1.74, (35) 

compared to the previous value of 1.91 [5] given in the Table I. After the publication of Temkin’s 
formula (34) and his numerical result in Eq. (35), Schwartz [7] did use such a slowly vanishing part 
1/r 2 in his trial wave function, obtaining precision results: 

a s = 5.965 ± 0.0003 and a t = 1.7686 ± 0.0002, (36) 

for the singlet (s) and triplet (t) scattering lengths. 


Symmetric Euler-Angle Decomposition 

In most of the scattering and eigenvalue problems, it is necessary to write the required wave 
function in terms of a product of two components: radial and angular. Only then the Hamiltonian 
operating on such a wave function gives equations which depend on the internal coordinates only, 
with terms arising from the operation of the angular part contained in the Hamiltonian. For 
one particle scattering from a fixed center, (or the relative motion of two particles) the Laplacian 
operator is of the form 


V 2 


— — ( — * a d 1 d 2 , 

r 2 dr V dr r 2 sin# dO U dO sin 2 # d(j) 2 


(37) 


where the radial and angular parts can be separated and the second term is propotional to the square 
of the angular momentum operator. When acting on a wave function which is an eigenfunction of 
angular momentum /, the Laplacian simplifies to 


v 2 = ——(r 2 —) - l SLLll 
r 2 dr ' dr r 2 ’ 

thus converting a 3-dimensional partial differential equation into an ordinary differential equation. 

In the case of two particles in an external force field, e.g., the field of the fixed nucleus (i.e., 
of infinite mass), the decomposition of the Laplacians is not as simple as in the case of relative 
motion of two particles. The wave function describing these particles is a function of the internal 
(i.e., meaning coordinates depending on the position of the nucleus and on two electrons in a plane) 
coordinates (ri, 1*2, #12) or (rq, 7*2, 7*12) and the angular component which is a function of three Euler 
angles which describe the orientation of the instantaneous plane formed by the two particles and the 
nucleus in space. The Euler angles are not unique. Breit [12] used the Euler angles which Hylleraas 
[13] introduced. These are two spherical angles of one the particles an’d the angle between the 
•7*1 — z plane and rq — r 2 plane, the internal coordinates being (7*1, t* 2 , #12)- Breit introduced these 
angles for P - wave functions and - because they are not symmetrically defined with respect to the 
two electrons - it is not easy to generalize this decomposition for all angular momenta. Temkin 
introduced a different set of symmetric Euler angles which allowed the separation of the radial part 
and angular part for any angular momentum /. The analysis of this problem was carried out by 
Bhatia and Temkin [14]. Fig. 1 contains a perspective drawing of the Euler angles which define 
the particle plane with respect to the space fixed x, y, and z axes. 
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The rotated axes x , y\ z with respect to the space fixed axes y, z are defined by 

f\ X r 2 


Having defined z\ define x 


Having defined z and x , define y 


The Euler angles are then 


z — 


|ri X r 2 I' 


x — 


| z x £ ; | 


j ^ / 

y — z X x . 


0 = angle between z and z , 


$ = angle between x and x, 


= angle between x and (r 2 — ?i). 

The operation of parity (ri — ri and r 2 -r 2 ) only involves the Euler angles 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


0 -+ 0 , 

S' 7T + $. (45) 
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Exchange corresponds to the transformation ri r 2 and is given by r\ V 2 plus the following 
transformation of the Euler angles 


0 -* 7T - 0, 

<J> -* 7T + <J>, 

'V -* 2tt - (46) 

The eigenfunctions of total angular momentum are purely functions of the Euler angles. Temkin 
has coined the name ’’rotational harmonics” for these eigenfunctions (they are usually called V 
functions). The angular momentum properties of these functions are 

M 2 p ; m - fc = ni(l + l)V™’ k , (47) 

and 

M z V™’ k - hmV™’ k . (48) 

They transform under parity as 

VV™’ k = V™' k (e,$, IT + V) = (-l)*p™’ fc (0,$,^), (49) 

and exchange as 

£uV™’ k {e, tf) = V™’ k (ir - 0, tt + $, 2tt - tf) = (-1 ) ; X>™’“ fc (0, $). (50) 

We see that whereas the V functions are eigenfunctions of parity, that they actually change 
indices under £ 12 . But operation of exchange commutes with the total Hamiltonian, therefore it 
must be possible to construct eigenfunctions of exchange: they are linear combinations of the above 
V functions. For m=0 they are 


+ _ p?* + (-1 
1 V5[i+w5-i)r 

and 


V 


AC — 

l 


V 0,K 


lV2 


Under exchange these linear combinations have the property 


£i 2 P ( k± = ±{-l) l+K V? ± . 


(51) 


(52) 


(53) 


Thus they are indeed eigenfunctions of exchange. Having constructed eigenfunctions (of the angular 
part) of exchange, one can construct a total wave function which has the correct properties under 
operations of parity and full exchange 


9 



Symposium on Atomic & Molecular Physics 


'Mri, r 2 ) = £[/f+(n,r 2 , 0 12 )V?+ (©, <1>, 4>) + //*" (n, r 2 , (0, 4>, )], (54) 

where the double prime on the summation indicates that the sum goes over every second value of 
k. The operation under exchange on Eq. (54) gives, 


fi 2 ®Jo(ri,r 2 ) = ® /0 (r 2 ,r 1 ) = £[/f+ (r 2 , r 1 ,0 12 )(-l) / +^f+(0, ¥) 

+ /r(»’2,ri,tfi 2 )(-l) ,+ " +1 ^r(©.^« r )]- (55) 

If therefore 

/r + (r2,ri,«i 2 ) = ±(-l) ,+ “/r + (ri,r 2 ,tfi 2 ), (56) 

and 

= ±(-l) ;+K+1 /r(ri,r 2 ,M. (57) 

then under complete exchange 

^i 2 ^/o(ri,r 2 ) = ^; 0 (r 2 ,ri) = ±®/ 0 (ri,r 2 ), (58) 

i.e., the total wave function will be symmetric (singlet) or antisymmetric (triplet). Under parity 
we see from Eq. (54) that 


Wio(r u r 2 ) = x f/o(-ri, -r 2 ) = (-l) K ^/ 0 (ri, r 2 ), 

so that parity is determined by evenness or oddness of k in Eq. (54). 

The kinetic energy part is given by 


(59) 


vi + v 2 = 


Id 2 Id 2 
ri dr\ Vl r 2 Or^ 2 


+ 


1 , 1 0 0 1 d 2 

r\ sin#! Mi "Mi sin 2 *! Oft’ 

d 1_ _1 . _0_ 1 d 2 

r 2 sin # 2 5 ^ 2 Sin 2 002 sin 2 0 2 dft 


(60) 


Historically, this was the extent of the development when I came to work with Dr. Ternkin 
in 1963. I worked for almost a year to convert the kinetic energy into Euler angles (plus internal 
coordinates) form. The following is the result 


vi + v 2 = 


1 o 2 


n + 


ri Oft r 2 Or 


0 2 . 1 
„2 r2 + ( „2 


+ 3 )( 


1 


0 


r 2 n sin^i 2 O0i2 


•sin# 


12 


0 

00 


12 


0+3- 


r 2 


(61) 
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Fi = 


^[sin 2 (9> + + cos ! (* + 5«.2)cot0^ 

+ cos ! (<t + 5*12)^555 + sin(2<f + 

1 cv 


sin \y v'y sin© d$ 

Si", 2 * + « 12 ) _i_ _|L + sin(2 * + ^JcotegAg 

2 1 . cot© d 2 , d 2 d 2 d 

2cos (* + 2^ 12 sin© + Al ^ + Bl #ip 

Ai = i + (cot 2 ©/sin 2 (?i 2) cos 2 (W + ^0i 2 ) 


( 62 ) 


(63) 


Bi — (cos$/sin 2 0i2)sin(^ + #12) - [sin(24> + 0i2)/sin 2 0i2sin 2 0] - -cot(-#i2). 


(64) 


The expressions for F 2 , A2, and B 2 can be obtained by replacing #12 by -#12 in the above expressions 
for t\ . Ai, and By . 

With these results the Schrodinger equation 


HVim = E^ tm 

for any m, can be reduced to radial equations which are independent of m: 


(65) 


[Lh 2 + 

+ 


2m f 


(. e - v W r - (4 + 4)[{ 'g±4_iir! + £ _ p«iL l{l + 

v ■’i Vf r ^i«. 2sin 2 6> 12 4 4sm6>i2 


fS^r«/r> + + ^ (1 _ *. _ *. + 

/l 1 \ r / 1 ^ 5 \ 1(1 ~h 1) f rAC _ 

+ - ^)[ _K (9 Cot ^ 12+ ^7 )f i “ 


12 


4sin^i2 


+ 


B ? +2 


/(«+ 2 )- 1 — _ 5 _ s _ s 2K )Bi K fi K ~ 2) ~] = 0 

4sme n Jl 4sin0 12 ' ° K lK 21 1 J 


(66) 


and 


2 m 


[l S i2 + -^{ E -v)]fr *)[{ 


1 1 uf /(/+ 1) - K 2 


2sin 2 ^i2 


ft 2 COt^i9 

H 1 — 

4 4sin#i2 


cot0 


+ tt-^(1 - 6 OK )B?+ 2 fi K+2) - + - 5 0k - 5 1k - 82 K )B lK f} K -*>~) 


12 


/(/ + l) 5 i K }f* 

(ac— 2)- 


4sin^i2 

+ L2 2 


4sin^i2 


1 




(1 - Son) 


4sin0 


12 


B? +2 t 


+ 2 )+ 


+ 


1 


-(1 — 5o« “ + ^2 nBlnfi ^ + ] — 0 


4sin0 


12 


(67) 
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where Lq u is the 5-wave part of the kinetic energy, and is the only term which survives in the 
description of S - waves 


and 




1 d 2 


r r 1 + 


ri dr 2 r 2 dr% 


d 2 ,1 1 , 

r 2 + (l2 + I 2 ) 


1 


5 


sin^io d0 


712 




12 


(68) 


Bln - 


[(/ — K + 1)(/ — K + 2)(Z + «)(/ + K — 1)] 2 
[1 + <5 2 *(\/2 — 1)] 


(69) 


Bf = B/ K [l + 5 2 «(\/2 - l)] 2 . (70) 

The above equations can be written in terms of (ri,r 2 ,ri 2 ) [14]. These equations have been 
generalized to the case when the nucleus is of finite mass [15], and also to the case when all the 
particles are of unequal masses [15]. 

The equation given by Breit [12] for P - wave can be obtained from our equations by noting the 
relation between his angles and our Euler angles (cf. Appendix of our paper [14]). 

I believe this analysis is one of the major mathematical achievements of Dr. Temkin and I am 
happy to have been a part of it. And this laid the mathematical foundation of much of our future 
research, also with Dr. Drachman. 


Optical Potential Approach for Scattering 

We follow here the Feshbach projection operator formalism [16] to obtain equations for the 
scattering function giving phase shifts in the elastic region. The phase shifts obtained have property 
of having a rigorous lower bounds. 

In this formalism, in order to project out the ground state we use symmetric projection operators 
P and Q, which for the hydrogenic (i.e., one-electron) target can be written expicitly [17] 

P = Pl+P2~PlP2 , (71) 


and 


Q = l-P, (72) 

which implies that P-\-Q — 1 and where the spatial projectors are such that for any arbitrary function 

$L(ri,r 2 ) 

Pi$L(ri,r 2 ) = ^o(ri) J <f> 0 {r , 1 )^ L (r , 1 ,r 2 )dr' 1 . (73) 

Note, Pi and P 2 commute with each other and are each idempotent, hence the complete P and Q 
operators are idempotent (P 2 =P; Q 2 —Q ) and orthogonal (PQ= 0). In the lim r — y oo, 


12 
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P^L -* sin(kr - + %), 

therefore P^l represents a scattering function, and 


(74) 


QVl -> 0. (75) 

Effectively Q^l describes the short range part of the total wave function of the system The 
Schrodinger equation is reduced, via an analysis which by now is well known [16], to an equation 
for P^l- 


[PUP + FH !? - E]P*l = 0 


(76) 


E-QHQ 

The middle term in Eq. (76) is the formal, but well defined expression for the optical potential V op . 
The total spatial function for the e-H and e-He + for the Lth partial wave is written as 


*L(ri,r 2 ) - ^^r LO (ri)^o(r 2 ) ±(lH 2) + * L (n, r 2 ), 


r 1 


(77) 


and 


P^ L = V L 

where the target function is given by 


(78) 


exp(-Zr 2 ). (79) 

The upper and lower signs correspond to singlet and triplet states, respectively. The first two 
terms containing ul explicitly give rise to the exchange approximation [1] and the function 4 >l is 
the correlation function. For arbitrary L this function is most efficiently written in terms of the 
symmetric Euler angles [14]: 4>£ has exactly the form of in Eq. (54), whereas here the open 
channel part is added explicitly in Eq. (77). 

The /’ s we here take as functions of rq, r 2 , and ri 2 . The ui(ri) of the scattered electron in Eq. 
(77) is determined by projecting on < Ylo(^i) ^ 0 ( 1 * 2 )’ 



/ 


Y L0 (il 1 )M^)(PHP + 


PHQ >< QHP 

E-QHQ 


- E)Pq L 


d£l\dv2 — 0 , 


(80) 


Carrying out the integration leads to an integro-differential equation for the scattering function 
ui(ri) and letting rq=r, 


dr 2 


L(L + 1) 
r 2 


+ v st (r) ± V ex - V op + k 2 


u L - 0 , 


(81) 


where v st (r) is the direct potential [Eq. (29)] and V ex is the non-local exchange potential of the 
“exchange approximation” [1] (cf. Eq. (28)). It should be noted that the many-body problem has 
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been reduced to one-body problem and the Eq. (81) can be solved for ul easily. The integral of 
the optical potential acting on t^(r) is: 


V op u L = r ( Y L0 <t> 0 PHQ — q HQ QHP9l) ■ 

The optical potental is expanded in terms of the eigenspectrum of the QHQ problem: 


(82) 


£ 


< > 
< > 


= 0 . 


(83) 


This leads to eigenfunctions and eigenvalues S s . By inserting a complete set of the functions 
obtained from the above equation into Eq. (82), V op ul can be written as 


Nw 

VopUi(ri) = 

s 


< YLoirJM^^Q^ >< > 

E-S s 


(84) 


TABLE II. e-H phase shifts of 1 5, 3 5, X P, and 3 P states for various k obtained from the method 
of polarized orbitals and from the optical potenial approach. 


Partial Wave 

k 

PO a 

OP 6 

PO a 

OP b 


Singlet 


Triplet 



0.1 

2.583 

2.55358 

2.945 

2.93853 


0.2 

2.144 

2.06678 

2.732 

2.71741 


0.3 

1.750 

1.69816 

2.519 

2.49975 

s 

0.4 

1.469 

1.41540 

2.320 

2.29408 


0.5 

1.251 

1.20094 

2.133 

2.10454 


0.6 


1.04083 


1.93272 


0.7 C 

0.947 

0.93111 

1.815 

1.77950 


0.8 C 

0.854 

0.88718 

1.682 

1.64379 


0.1 

0.0067 

0.0063083 

0.0109 

0.010382 


0.2 

0.0171 

0.014988 

0.0486 

0.045345 


0.3 

0.0210 

0.016613 

0.1151 

0.10679 

p 

0.4 

0.0163 

0.099980 

0.2005 

0.18730 


0.5 

0.0064 

-0.00084017 

0.2867 

0.27058 


0.6 

-0.0039 

-0.010359 

0.3574 

0.34128 


0.7 

-0.0100 

-0.013483 

0.4063 

0.39257 


0.8 

-0.0095 

-0.0048524 

0.4351 

0.42730 


a S - wave phase shifts obtained by Temkin and Lamkin [5] using the method of polarized orbital. 
The P-wave results are from Sloan [6]. 

b Phase shifts obtained from the optical potential approach. S - wave results are from Ref. [19] and 

P- wave results are from Ref. [20]. 

c The polarized orbital results are from [9]. 
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For 5-wave (i.e., L=0), X>£=constant and the correlation function is only a function of the radial 
coordinates. The formalism up to this point had already been developed. The innovation which 
Dr. Temkin and I introduced was the use of a correlated function of Hylleraas form 

d > L=0 = e -7n £ C lmn r[r?r? 2 ± (1 # 2 ) . (85) 

linn 

In particular the projection of Pi on 4>£o requires a nontrivial analytical integration [18]. Here the 

sum includes all triples such that l + m + n=u> and u= 0, 1, 2, , 7,8,9. The number of terms for 

each uj depends on spin and whether y—8 or not. For P- wave (i.e., L— 1) the correlation functions 
are also taken of the Hylleraas form. 

To summarize the calculations, the QHQ problem is solved (for a given L, 7 and 5 and N u ). 

The result is a set of eigenvalues £ s (s= 1,2, NJ) and associated eigenfunctions From 

them the optical potential, Eq. (82) is constructed, and the integro-differential Eq. (81) is solved 
noniteratively. The solution is unique (up to an arbitrary normalization) with asymptotic form 

lim UL(r) oc sin(kr - + Z 1 ln( 2 kr) + arg[r(l - ^ )] + rji). ( 86 ) 

Electron-hydrogen phase shifts (i.e., Z=l) are given in Table II. For 1 5-wave scattering, the 
optical potential (OP) phase shifts converged to the accuracy shown when the maximum number 
of terms in the wave function was 95 and in 3 5 -wave the convergence was obtained when the 
number of terms was 84. The P- wave phase shifts converge slowly compared to the 5-wave results. 
Therefore, the computation was carried up to 220 terms. The convergence [19, 20] of results suggest 
that they are accurate up to five significant figures after the decimal and to that accuracy they 
are rigorous lower bounds, provided the total energy of the system is less than those of all the 
resonance positions [21] . Phase shifts are compared in Table II to the polarized orbital results 
of Temkin and Lamkin [5] for 5- wave, for P- wave with those of Sloan [ 6 ], including the exchange 
polarization terms. The polarized orbital method does not provide any bound on the phase shifts 
but they are seen to contain the dominant part of the correlation enhancement over the exchange 
approximation. In particular the polarized orbital l P results show the correct undulations (as a 
function of k ) as the precision results. 

A similar calculation has been carried out for the scattering of electrons from helium ions [22,23]. 
The non-Coulomb part of the phase shifts as a function of k are given in Tables III for l S , 3 5, l P 
and 3 P. In this case because the Coulomb field extends very far, Eq. (81) has to be integrated to 
large distance especially for small values of k. 

In the singlet P case, the exchange approximation results are negative and there is a lot of 
cancellation with the contributions to the phase shifts from the optical potential, unlike in the 
triplet case where the exchange approximation results are always positive. 
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TABLE III. e-He + phase shifts of 1 5, 3 5, X P, and 3 P states for various k obtained from the 
method of polarized orbitals and from the optical potenial approach. 


Partial Wave 

k 

Exch 

PO“ 

OP 6 

Exch. 

PO a 

OP 6 



Singlet 



Triplet 



0.4 


0.4301 

0.42601 


0.9235 

0.91300 


0.5 



0.41964 



0.90275 


0.6 


0.4153 

0.41278 


0.9015 

0.89050 

5 

0.7 

0.34871 


0.40561 



0.87640 


0.8 

0.33949 

0.3986 

0.39857 

0.85195 

0.8723 

0.86069 


0.9 

0.33100 


0.39202 

0.83444 


0.84356 


1.0 

0.32302 

0.3823 

0.38634 

0.81645 

0.8371 

0.82531 


1.1 

0.31577 


0.38187 

0.79764 


0.80625 


1.2 

0.30940 

0.3685 

0.37899 

0.77821 

0.7984 

0.78666 


1.3 

0.30400 


0.37832 

0.75851 


0.76684 


1.4 

0.29957 

0.3579 

0.38560 

0.73872 

0.7591 

0.74697 


0.1 

-0.070322 


-0.038311 

0.17705 


0.21516 


0.2 

-0.074161 

-0.0404 

-0.038958 

0.17913 

0.2232 

0.21681 


0.3 

-0.075219 


-0.039911 

0.18238 


0.21944 

P 

0.4 

-0.076473 

-0.0428 

-0.040971 

0.18653 

0.2290 

0.22283 


0.5 

-0.077742 


-0.041951 

0.19123 


0.22661 


0.6 

-0.078820 

-0.0450 

-0.042633 

0.19615 

0.2364 

0.23048 


0.7 

-0.079501 


-0.042834 

0.20096 


0.23415 


0.8 

-0.079626 

-0.0457 

-0.042383 

0.20540 

0.2429 

0.23744 


0.9 

-0.079078 


-0.041158 

0.20927 


0.24008 


1.0 

-0.077762 

-0.0436 

-0.039036 

0.21245 

0.2469 

0.24202 


1.1 

-0.075639 


-0.035948 

0.21489 


0.24322 


1.2 

-0.072707 

-0.0384 

-0.031838 

0.21659 

0.2479 

0.24378 


1.3 

-0.068996 


-0.026592 

0.21757 


0.24370 


1.4 

-0.064574 

-0.0301 

-0.019982 

0.21790 

0.2317 

0.24320 


“S'-wave and P-wave phase shifts obtained by Sloan [6] using the method of polarized orbital. 

6 Phase shifts obtained from the optical potential approach: S - wave results are from Ref. [22] and 
P- wave results are from Ref. [23]. 


Resonances in Two-Electron Systems 

In the above section, projection operators P and Q for one-electron targets were given. These 
have been used to calculate resonance parameters for a number of 1,3 5, 1,3 P, and 1,3 P states in H“ 
and He. Unlike other methods where one has to hunt for resonance positions, they are obtained by 
optimizing the functional < &QHQ& >/< > by using the Rayleigh-Ritz variation principle. 

The positions obtained do not include the shift due to their being embedded in the continuum [17] 
and this shift is calculated separately using various approximations for the continuum functions. 
For illustration, only two sets of results (positions and widths) are given in Table IV for X S in H” 
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Table IV. Some results for resonance states in H and He 


State 

Position (eV) 

Width (eV) 

Position (eV) 

Width (eV) 

^(H-) 

9.55735±0.00005 a 

0.0004717±0.00002 

9.558±0.010 6 


1 P(He) 

60.133±0.015 c 

0.038±0.004 

eo.mio.oio'* 

0.038±0.002 <i 


a Resonance position is with respect to the ground state of H atom [24] 
Experimental result are from Sanche and Burrow [27]. 
c Resonance position is with respect to the ground state of He atom [25]. 
Experimental results are from Morgan and Ederer [28]. 


[24] and X P states in He [25]. This odd-parity resonance is the lowest one in the series of resonances 
observed in vacuum ultraviolet absorption by Madden and Codling [26]. 

We see that results have been calculated with high precision and they agree with the experimental 
results. The line shape parameter ^=-2.80±0.025 which we have precisely defined and calculated 
by further manipulation of the Feshbach theory [25] also agrees very well with the experimental 
result -2.55±0.16 [28]. 


Projection Operators for More Than Two-electron Systems 

For targets with more than one electron, it is difficult to construct projection operators which 
are idempotent, i.e., P 2 =P and Q 2 —Q. One of the difficulties is that target wave functions for 
more than one-electron targets cannot be written down exactly. But we can construct them in such 
a way that the matrix elements < $ rn Q 2 $ n > =< ^ m Q^ n > for any arbitrary antisymmetric wave 
functions and 4> n . One of our motivations was to study 2 S resonance below the 3 5 (elastic 
region) of the helium atom. 

In order to include all coordinates explicitly, we assume LS (i.e, Russel-Saunders) coupling 
and introduce channel functions in such a way that the ground state <fto is coupled to the angular 
momentum and spin \ of partial wave of the incoming electron 

Mr (t) ) = ^(L 0 / ! Mom i |LM)(5 0 tMs 0 m s |5M s )^o(a ; W)y / , ro ,(a)Xi ros (*)- (87) 

In Eq. (87) (both space and spin) indicates the absence of the ith coordinates from the total 
(N+ 1) coordinates in the electron-target system. The a;M signifies 

^ — (*E«-j-l ? — i) . (88) 

Eq. (88) implies that the target has N electrons and rW in Eq. (87) implies the absence of r 2 - 

rW = (fii,^;x (,,) ). (89) 

Let pi represent a cyclic permutation, so that 
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(— l) Pi — parity of the permutation which carries [1, 2, ..N + 1] into [i, i + 1, ..i — 1], (90) 

Specifically (-l) Pt =l for (iV + l)=odd, (-l) p *=(— 1) 2_1 for (JV + l)=even, and (—l) Pl =l for all N. 
Following Feshbach [16], we can define the P^i part of the total wave function as 

N+l 

( 91 ) 

i= 1 

Here u\(ri) are scattering functions which have not been specified but have the asymptotic property 
for ri -» oo 


lim ri _> 00 ui(r i ) = 


sin(krj — hr/2 + r]\) 

kn 


(92) 


which implies that both P'S* and ^ have the same asymptotic form 


llm n ^P« = = ( _ 1)Pi sin(kn-l,/2 + ^) ^ (rW) 


(93) 


To derive a specific form of P (and Q) we also require that Q ^ have no ground state in it for 
any coordinate of the scattered ri (not only as r t — oo) 


< fo(r®)(l-P)9 >= 0, 

which for the purpose of the derivation can be written 


< '0o (r^)® >=< Mr (t) )P* > • 


Define the left hand side as 


w(ri) = (-l) p * < i t/> 0 (r w )W > r w • 
Substituting Eq. (91) into Eq. (95), we can express Eq. (96) as 


/*CO 

w(r,) = u(ri) - / K(ri\rj)u(rj)r 2 jdrj 
J 0 


where the kernel K is given by 


(94) 


(95) 


(96) 


(97) 


K(n\ rj ) = (-1 )Pi+Pj+'N < Mr (i) )Mr (j) ) > r W) ■ (98) 

Here the integration is over which denotes all coordinates except r{ and r 3 (i ^ j). The kernel 
K can be expressed in terms of discrete and orthonormal set of eigenfunctions v a 


V a {ri) = A a f 

Jo 


K(r i \r2)v a (r 2 )rldr 2 . 


(99) 
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Some reflection will show that K(ri\r 2 ) can be written as 

K (r,|r 2 ) = £ (100) 

then substituting Eq. (100) in Eq. (99) the latter will be an identity. Using Eq. (100) in Eq. (97) 
allows u(r z ) to be written 

«(r.) = £ ».(ri) < »»« > +»(ri) + £ < T > ■ <U>1) 

xS. A % X D ~ 1 

Here vp is the eigenfunction associated with A p and v a with A a =l. It can be verified that the first 
term in Eq. (101) does not contribute to the projection < > (cf. Refs. [29,30]). Thus we 

can write 


< <Mr (1) )P» >=< 4 11 £ (-irMr.) + E > ]0o(r< i >) >,<„ . (102) 


A/3 


Substituting for w(ri) from Eq. (97) and rearranging, we get 


<*o(r (1) )Ptf >=< rf’l ><tf'i)+E ,,( 1 < r- > ]»>' nos) 

i— 1 A a 


Aa-1 

By comparing the left- and right-sides of Eq. (103), The expression for P can be extracted 


P = X>o(r {!) >< Mr (i) ) + E Va{ri)Mri l ) l < (104) 

.• = 1 Aq ^ 1 


which implies 


r\ 1 r / / i V X /i n r\ 

Q = 1 - 2^ m >< V>o + 2 ^ — t i — ]• ( 105 ) 

i= 1 Aa Aq! 1 

These operators are symmetric in all (N + 1)- -particle coordinates and they contain the dependence 
on space and spin of all coordinates explicitly. The main advantage of these expressions for P 
and Q, as compared to the heuristic form given by Feshbach, are they are complete and explicit. 
In addition they have been extended to scattering (including resonances) occuring in the inelastic 
domain (cf. Appendix B of Ref . [29]). 

It has been shown in Ref. [30] that even though the above operators are not idempotent 
(P 2 ^ P and Q 2 ^ Q) as operator identities that the matrix elements are equivalent < $;P 2 <I>j > = 
< $iP$j > and < >— < > for any arbitrary antisymmetric functions <I>; and $j. 

When the second term in Eq. (104) and the third term in Eq. (105) are dropped, we get 
the quasi-projection operators P and Q. These operators were employed by Temkin et al. [31] 
to calculate the resonance parameters of the He“ [ls(2s) 2 ] 2 S resonance (first observed by Schulz 
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Figure 2: Metastable helium excitation from Schultz and Fox [35]. 


[32]), using open shell and closed shell target wave functions <j>o and an angle-independent trial wave 
function as well as a configuration-interaction wave function containing up to 40 configurations. 
They obtained E res ( 2 S) =19. 363 eV plus a width r=0.014 eV. The difference between the results 
obtained by using the two different wave functions </>o is of the order of 0.02 eV. The resonance 
position agrees with the experimental result 19.31±0.03 eV of Kuyatt et al. [33]. Calculations have 
also been carried out by Berk et al. [34] using target wave functions going up to a 10 term Hylleraas 
type wave functions and the full projection operators P and Q given in Eqs. (104) and (105). Again 
configuration interaction type wave functions containing up to 40 terms were employed obtaining 
resonance position of (1.s2s 2 ) 2 S He - which is ~0.013 eV above the experimental position at 19.37 
eV given by Schulz and Fox [35]. (That difference is presumably due to the uncalculated shift, 
which also occurs in the Feshbach theory.) 

Now we come to 2 P wide resonance in He - above the 2 3 5 threshold but below the 2 l S threshold 
of He, which has been first observed by Schulz and Fox [35] (cf. Fig. 2). 

Initially, it was thought to be a shape resonance because of its being above the 3 S threshold. 
Bhatia and Temkin [36] used the same program which was used to calculate the 2 S resonance below 
the 3 S threshold and where quasi-projection operators were employed. Only those configurations 
in the trial wave functions were included which were orthogonal to the 3 S state of He, thus avoiding 
the need to project out the 3 S state. The calculations were carried out by using closed shell as 
well as open shell functions in the projection operators, giving the resonance position at 20.52489 
and 20.56029 eV, repectively, for a trial wave function consisting of 40 terms. The position agrees 
with the result 20.536 eV obtained by Chung using his hole-projection technique [37]. The partial 
widths to 1 S and 3 S thresholds were found to be 0.0024 and 0.437 eV, and the calculation also 
showed that the total width is dominated by the decay to the (excited) 2 3 S state of He. These 
results also indicated that the resonance was a Feshbach resonance associated with the closed 2 l S 
state of the target He rather than a shape resonance caused by the open 2 3 S state. An accurate 
calculation carried out by Junker [38] using the complex rotation method gave 20.33 and 0.575 
eV for the position and width, respectively. These results agree with average experimental results 
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20.3dh0.3 and total width 0.5 eV of Schultz and Fox [35], and Brunt et al. [39]. 

I shall not discuss here Dr. Temkin’s (and collaborators’) work on dispersion relations [40, 
41] and in particular his work on threshold laws for electron-impact ionization of atoms and ions 
[which also apply to photon double (detachment/ionization) of (negative ions/atoms)]. The latter 
is discussed by Dr. R. Wehlitz in these Proceedings. A dispersion relation (DR) relates the real part 
of the scattering amplitude to an integral over the imaginary part, which in turn is proportional to 
total cross section. Such relations are important in judging the consistencies and accuracy of both 
theoretical calculations and experimental measurements. The problem in electron-atom scattering 
has been to correctly include the effects of exchange. Dr. Temkin has proposed using partial wave 
dispersion relations to solve this problem. At this point the correct partial wave DR have been 
constructed both approximately [42] and exactly [43] in the static-exchange approximation. 

Up to now, I have described some of Dr. Temkin’s important contributions in electron-atom 
scattering and associated problems in atomic physics mentioned above . He has also made impor- 
tant contributions in the field of electron-molecule scattering: Briefly stated, Dr. Temkin (with 
various coworkers) introduced the ’’fixed nuclei” approximation (as well as the name) in elect ron- 
(diatomic) molecule scattering [44,45]. They showed that the (partial wave) scattering amplitude 
could be expressed as the product of the two factors; one of which depends on scattering parameters 
resulting from the dynamical interaction of the electron with the molecule, which is most conve- 
niently calculated in the body-frame of the molecule, and a second factor, depending on geometrical 
functions representing the rotation of the scattering angles from the body to the lab frame. The 
cross sections, averaged over orientations of the internuclear axis could then be expressed as a sum 
of scattering parameters multiplied by spherical harmonics together with vector coupling coeffi- 
cients. A second contribution concerns the ’’adiabatic nuclei” approximation, originally introduced 
by Chase [46] in the context of nuclear physics. Chase showed that in a cogent approximation that 
the amplitude for rotational excitation can be expressed as a matrix element of the fixed-nuclei 
amplitude between initial and final rotational states. Because the dependence on the orientation is 
analytic (actually V function), these integrals can also be done analytically, again - first - by Temkin 
and coworkers [46,47]. (The name adiabatic-nuclei was also coined by Temkin [48].) Temkin was 
also involved in other developments, most notably the hybrid theory [49], which will further be 
discussed by Dr. B. I. Schneider in his contribution to these Proceedings. 


Having described some of Dr. Temkin’s work, I now will describe some of the work carried out 
by Dr. Drachman. Perhaps, after Sir Harrie Massey, Dr. Drachman has not only made important 
contributions to positron physics but has also made it a popular subject of research. 


Scattering of Positrons from Hydrogen Atoms 

Calculation of positron-hydrogen scattering at low energies can be carried out by the method of 
polarized orbitals as in the case of e-H scattering. Instead, Drachman [50], employing a variation 
of the method, chose the wave function of the form 


tf(ri,r 2 ) = w(r 1 )[l + G(ri,r 2 )]^ 0 (r2). 
The function G(ri, r 2 ), correct to first order in the potential 


(106) 
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V(ri,r 2 ) = t 7 , 

r 1 |r 1 -r 2 | 

has been determined by Dalgarno and Lynn [51] as a solution of the differential equation 


(107) 


[G, H (2)]</> 0 (r 2 ) = (V- < V >)^o(r 2 ) (108) 

where 

H{2)<j> 0 {r 2 ) = [-V 2 - -fo>(r 2 ) = -0o(r 2 ), (109) 

r 2 

and 

< V >= J dr 2 (j>o(r 2 )V ( r i) r 2 )^’o( r 2 )- (110) 

The function G includes all multipoles >0. The adiabatic equation, correct to second order in V , 
is given by 


[— vf - k 2 + + F 2 ]u(r x ) = 0. 


(HI) 


In the derivation of the above equation, we have taken < G >=0. Here 

Vst = 2e- 2r '(l+-), 

r l 


( 112 ) 


and 


V 2 = rf 2 [ 5 - (4r? + 8ri + 10)e- 2ri + (4r? + 7r? + 8r x + 5)e~ 4ri 
-2(r x + l) 2 (e -2ri + e -4ri )(Ei[2ri] - 21 n[ 27 r x ]) 

— 2Ei[— 2r 1 ]([r 1 - l] 2 e 2ri + [r 2 + 2r x - 3] + 4[n + l]e- 2r ‘)], (113) 

where ln7=0. 57721 is the Euler’s constant and 

r oo p -y 

Ei[— z] - - / — dy. (114) 

j z y 

For small value of ri, v si 2rj” 1 -2 and V 2 — > -1, while for rq — > 00 V2= -4.5r“ 4 . Eq. (Ill) 
can be solved for various values of k to determine the phase shifts and the results are given in 
Fig. 3. For k= 0, the scattering length requires a correction due to the long-range potential [11], as 
indicated in Eq. (34), and its value is -2.54 versus Schwartz’s -2.10 [7], which is the upper bound 
on the scattering length. This shows that the potential in Eq. (113) is too attractive. Drachman 
[50] modified the potential to 


v s t + V 2 + (a — l)V 2 o, 


(115) 
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Figure 3: Phase shifts in radians for p- and d - wave scattering. The solid curves are for a?=0.1, and 
the vertical bars indicate a=1.0 and a=0 limits. The triangles are Brandsden’s p - wave values (Ref. 
[52]). 


where V 20 is the monopole part of V 2 . The complete suppression of monopole or short-range part 
of the potential, i.e., a=0 gives scattering length -2.07, in good agreement with Schwartz’s value. 
Drachman calculated P - wave and D - wave phase shift for various values of k and found reasonable 
agreement with those of Brandsen [52] (cf. Fig. 3). 

In the same calculation, Drachman concluded, from the change of sign of the scattering length 
with increasing mass of the positron, that there is no bound state of the system e+-e“-P unless 
the mass of the positron is 3.6 m e for 01 = 0 . 1 , while for a=l it is about 3.1m e . 

Houston and Drachman [53] using a more flexible wave function than that of Schwartz [7] in the 
Kohn variational method, obtained an upper-bound scattering length a < -2.10278. They obtained 
an extrapolated estimate ‘a— - 2.1036±0.0004. These results are in good agreement with Schwartz’s 
result a < -2.10. They applied the Harris method [54] to obtain 5-wave phase shifts at nonzero 
energies which are in good agreement with those obtained by the Feshbach formalism [16] described 
below. 

One of my first calculations with Drachman [55] was the 5-wave elastic scattering of positrons 
(e + ) from hydrogen atoms below the positronium pickup threshold. As mentioned above for 
electron-hydrogen scattering, rigorous lower bounds have been obtained using the Feshbach projec- 
tion formalism [16]. Similarly, the results for e + -H have rigorous lower bounds. Since there is no 
exchange between a positron (labeled 1) and an electron (labeled 2) in the hydrogen atom in this 
process, the projection operators P and Q are defined as 


P = |^o(r 2 ) >< 0o(r 2 )|, 


(116) 
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Q = i - P, 


(117) 


such that P operating on the wave function does not change its asymptotic form in the limit ri — > oo 




(118) 


TABLE V. e + -H phase shifts for 5-wave scattering for various k. 




S - wave 



P- wave 


k 

Schwartz [7] 

^extrap 

A T) 

^final 

Armstead® [8] 

^final 

0.1 

0.151 

0.148085 

0.000223 

0.1483 

0.008,0.009(1) 

0.0094 

0.2 

0.188 

0.187496 

0.000200 

0.1877 

0.032,0.033(1) 

0.0338 

0.3 

0.168 

0.167407 

0.000306 

0.1677 

0.064,0.065(1) 

0.0665 

0.4 

0.120 

0.119724 

0.000420 

0.1201 

0.099,0.102(1) 

0.1016 

0.5 

0.062 

0.061934 

0.000429 

0.0624 

0.130,0.132(1) 

0.1309 

0.6 

0.007 

0.003191 

0.000689 

0.0039 

0.153,0.156(2) 

0.1547 

0.7 

-0.054 

-0.052183 

0.000980 

-0.0512 

0.175,0.178(3) 

0.1799 


a Armstead [8] has given two sets of results. The first entry gives his converged results while 
the second entry his estimate of most probable results with uncertainty in the last figure given in 
perenthesis. 


Here T' is given by 


®l( ri,r 2 ) = UL(ri)<Mr 2 ) + $(ri,r 2 ), (119) 

where the generalized Hylleraas function is 

N 

* = e -(Sr l+ ^+ar n ) £ Clmn r[r?r? 2 . (120) 

Imn 

In the absence of exchange, the e+-H problem should be easy to solve but not so due to the 
virtual positronium formation. Therefore, we have included e ~ arL 2 in $, where ri 2 is the distance 
between the positron and the electron. An integro-differential equation of the form Eq. (81) 
{V e x = 0 here) is solved for the scattering functions ui(ri) and phase shifts t]l are obtained in the 
limit ri — ^ oo from 


UL{r i) = sin(kri - Ln/2 + 1 %). (121) 

It should be pointed out that the phase shifts are negative in the absence of the optical potential. 
The inclusion of the optical potential, which is attractive as in the e-H scattering, increases the 
phase shifts from the values obtained in the presence of only the repulsive static potential v s t in 
Eq. (81). The phase shifts were calculated for up to N=84 and extrapolated for N ^ oo. These 
phase shifts plus a correction At] for the long-range polarization potential are compared in Table 
V with those obtained by Schwartz [7] using the Kohn variational principle. The long-range effects 
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are not well included in a Hylleraas type correlation functions and therefore have to be calculated 
separately. Our final results (//final = ^extrapolated + A 77 ), which, we believe, are accurate within 
±0.0002 radians, differ at k= 0.6 and 0.7 from those obtained by Schwartz [7]. 

A similar calculation [56] has been carried out for P - wave scattering. Here the closed channel 
function is given by 


$ = /icos(^ 12 )P 1 1+ + f 2 sin(^ 12 )Pj-. (122) 

The fi are taken of the Hylleraas type with two non-linear parameters: 

/I=e -(V 1+ „) EEE 

^ Imn r l r 2 r l2i (123) 

/> 1 m > 0 >0 

and 

/2 = e -(w^ 2 ) EEE 

C\mn r 2 r?r? 2 . (124) 

/>1 771 >0 >0 

The fi are linearly independent functions because the positron and the electron are distinguishable. 
Now there are four nonlinear paremeters to be veried to get the best results. The maximum number 
of terms for each / is 84, giving a total of 168 terms in 4>. Here we have added the dipole adiabatic 
and the quadrupole plus nonadiabatic corrections to the extrapolated phase shifts. Our final results 
are compared to those of Armstead [8] in Table V. The rigorous lower bound is lost due to the 
addition of corrections for the long-range potential. Nevertheless, the results still are accurate. 

Needless to say, these results have stood the test of the time and are still considered to be the 
benchmark results. 


Annihilation of Positrons with Electrons, Z e ff 

An important process is the annihilation in flight of positrons by atomic hydrogen resulting in 
the 511 KeV line which has been observed in solar flares [57] and from the galactic center [58]. 
This line can be used to infer properties of flares and the solar plasmas. Having calculated the 
wave functions for the scattering of positrons from hydrogen atoms, the partial cross sections for 
annihilation can be calculated from the expression given below [59] 

<?a = Z ea ^(wa 2 0 ), (125) 

where a is the fine-structure constant, a 0 is the Bohr radius, and k is the incident positron momen- 
tum in units of aQ 1 . The quantity Z e //, whch is the measure of the probability that the positron 
and electron are at the same point, depends on specific properties of the positron-atom system. 
Z e ff is equal to Z, the number of electrons in the atom, when the positron can be represented as a 
free particle. For the hydrogen atom 


Z e f{- J j dridr 2 |'$’(ri,r 2 )| 2 6(ri — r 2 ), 


(126) 


25 



Symposium on Atomic & Molecular Physics 


where is the positron-hydrogen wave function. It should be noted that for the calculation 

of Z e //, ul(t) should have a plane wave normalization, i.e., 


« L (r) = (2 L + l) 3i "^ - r -- | r L ^ > PL{coaW) (127) 

We present Z e // for L — 0 and 1 [60,61] in Table VI and compare them with those obtained by 
Humberston and Wallace [62] and Humberston [63]. For L > 1, we use the plane wave expansion 
for the incident positron 


00 POO 

Z e f[{L > 1) = £(2L + !) J o drr2 ^o( r ) jl( kr ) ■ 

Using the identity 

oo 

£( 2L + l)j 2 L = l 

L = 0 


we get [with 4>o(r) = 2e r ] 


(128) 


(129) 


Z e f[(L > 1) 


POO 

4/ drr 2 e-^(l-jl-Zjl) 
J 0 


k 2 

k 2 + 1 



(^ ln (! + k2 ) 


1 + k 2 /2 
1 + k 2 ’ 


We give the total Z e r in Table VI. 


(130) 


Table VI. Z e j j for e + -H annihilation. k = 0 results are from Ref. [64,65]. 



O 

II 

h-} 


L= 

=1 

Z e ff (Total) 

k 

Z eff (L = 0) 

Ref. [62] 

Z e ff{L = 1) 

Ref. [63] 


0 





8.868 

0.1 

7.363 

7.5 

0.13008 

0.1335 

7.493 

0.2 

5.538 

5.7 

0.53994 

0.5366 

6.079 

0.3 

4.184 

4.3 

1.12441 

1.114 

5.312 

0.4 

3.327 

3.3 

1.76292 

1.719 

5.100 

0.5 

2.730 

2.7 

2.33910 

2.353 

5.091 

0.6 

2.279 

2.3 

2.84988 

2.823 

5.168 

0.7 

1.950 


3.67030 

3.637 

5.683 


We can calculate the thermally averaged annihilation parameter 

POO 

Z(T) = / dkf T (k)Z G ff(k), (131) 

Jo 

where fr{k) is the Maxwell-Boltzman distribution function. We can fit the calculated Z^(k) to a 
sixth-degree polynomial of the form 
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Figure 4: Thermally averaged annihilation parameter Z( T) obtained from the polynomial fit 
Zeff(k). 


6 

Z ef[ (k) = '£Z n k n . (132) 

0 

Using this expansion, the integration in Eq. (131) can be carried out analytically and we have 

Z(T) = Y,Zn(T4) n/2 , (133) 

71—0 

where T 4 is the temperature in units of 10 4 K and 

= a n Z n A~ n/2 . (134) 

Here Ao=15.789 and Z n , Z n , and a n are given in Table VII. 

Table VII. Coeffients of the polynomial, Z n , Z n , and a n . 


n 

z n 

Z n 

Qn 

0 

8.868 

8.868 

1 

1 

-7.38 

-2.226 

27T- 1 / 2 

2 

-102.77 

-9.763 

3 

2 

3 

527.38 

18.971 

47T- 1 / 2 

4 

-978.68 

-14.722 

16 

4 

5 

773.15 

5.284 

127T -1 / 2 

6 

-197.17 

-0.658 

105 

4 


Now Z(T) can be calculated at any temperature (cf. Fig. 4). 
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e"*~-He Scattering 

Although accurate calculations for e+-H can be carried out because wave functions of hydrogen 
atoms are known exactly it is difficult to perform experiments on this system. On the other hand, 
e"*~-He experiments can be performed relatively easily but the calculations are rather tedious because 
elaborate He functions are difficult to employ. However, Houston and Drachman [53] used those 
simple wave functions which gave a reasonable value of the polarizabilty of the He atom. They used 
the Harris method [54] to calculate phase shifts at low energies and their results agreed with the 
results of the variational calculation of Drachman [65]. They added a term to the wave function to 
represent the long-range dipole potential and obtained scattering length a=-0.524 and Z e f /— 4.3 at 
k 2 — 0, using the Kohn variational method. 


Properties of Ps and Photodetachment 

The positronium negative ion (Ps - ), consisting of two electrons and a positron, is particle stable 
and decays only by e + -e - annihilation into gamma rays. Mills [66,67] has produced and detected 
this ion and measured its lifetime. Drachman and I [ 68 ] calculated its ground state ( 1 5 f ) energy by 
using a trial function of the Hylleraas form, calculated expectation values of delta functions, and 
cusp conditions given by 


|C£) 

II 

< <5(r*) > \ 

(135) 


) < *(ri 2 ) > _1 • 

(136) 


Table VIII. Binding energy (Ry) of 1 S state of Ps , expectation values of S functions, cusp 
conditions, and decay rate T((nsec) -1 ). [The notation A(— B) stands for A X 10 -jB .] 


N 

7 

<5 

Binding energy 


S( r 12 ) 

v% 

v 12 

r 

125 

0.2585 

0.3585 

0.024009788 

0.020722 

1.7151(-4) 

-0.49910 

0.49711 

2.0850 

161 

0.3700 

0.3700 

0.024010026 

0.020732 

1.7136(-4) 

-0.49986 

0.49695 

2.0860 

202 

0.380 

0.380 

0.024010089 

0.0200730 

1.71 29 (-4) 

-0.49740 

0.49740 

2.0858 

120 

0.604 

0.296 

0.024009966 

0.020733 

1.7190(-4) 

-0.50000 

0.49347 

2.0861 

165 

0.604 

0.314 

0.024010079 

0.020773 

1.7164 (-4) 

-0.49999 

0.49441 

2.0861 

220 

0.604 

0.313 

0.024010113 

0.020733 

1.7150(-4) 

-0.50000 

0.49508 

2.0861 


Here ri and r 2 are the relative distances of electrons 1 and 2 with respect to the positron, and 
r 1 2 = |ri — r 2 |. Results are in given in Table VIII for the Hylleraas wave functions (see Eq. (141) 
below) with up to V=203 with 7 — 5, and up to 220 terms with 7 / 5, respectively. These 
functions have been used to calculate other properties. The cusp conditions test the accuracy of 
the wave functions near points of coalescence, since jq = i/ 2 = — \ and jq 2 = for exact solutions 
of the Schrodinger equation. The convergence of results, given in Table VIII, shows that our wave 
functions should be fairly accurate. To a sufficient accuracy the Ps - decay rate is given by 
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T = 2na 4 (c/ao)[l — a(17/n — 197r/12)] < $(ri) > 

= 100.6174 < 8(ri) > (nsec) -1 , (137) 

where the correction term proportional to a is due to the triplet lifetime [69] and the leading 
radiative correction to the singlet lifetime [70]. Our calculated value is in agreement with the 
measured [67] value T = 2.09 ± 0.09(nsec) -1 . 

Mills [71] recognized that the existence of 3 P e state of positronium, as in H - [72], would have 
interesting experimental consequences: The state would be metastable against breakup because 
Ps ~( 3 P e ) -> Ps^S) + e~ is nonrelativistically forbidden. Using 70-term Hylleraas wave function, 
Mills [71] did not find such a state . We [68] too failed to find this state even when we used larger 
expansions and also included long-range terms of several types [73]. 


Photodetachment of Ps 

Mills [66] suggested that Ps” ion could be used to generate positronium (Ps) beams of controlled 
energy; this would involve acceleration of Ps - ions and photodetachment of one electron. Drachman 
and I [74] calculated the dipole transition matrix elements by two simplifications: the intial Ps” 
wave function is represented by an asymptotic form whose normalization comes from our most 
accurate wave function [68], and the final state is a plane wave. 

The Hamiltonian of the system consisting of two electrons (pi, p 2 ) and one positron (x) is 


H = — - Vi + V* - 


p 2 


+ 


|/>1-X| |p 2 -x| \Pl-P2\' 


The center-of-mass system is 


R= g(Pi + P2 + x), ri=pi — x, r 2 = p 2 -x, 


(138) 


(139) 


where R is the coordinate of the center of mass of the entire system, rj. and r 2 are the distances of 
electrons 1 and 2 from the positron, respectively. The Hamiltonian is given by 



vi + vi + v ri 


• V 


r 2 


1 

+ — + 
r i 


»'2 


»*12. * 


(140) 


Omitting the center of mass coordinate R, which describes uniform motion of the center of mass, 
we write a wave function for the Ps” ground state in the Hylleraas form 


^(n,r 2 ) = *<(»•!, r 2 ,r 12 ) = -2= + r£r?*e-™-* r ‘K 2 . (141) 

V8^ 2 Imn 

The final state consists of an electron in a p state moving away relative to the center-of-mass of the 
Ps atom. We use the coordinate R 2 = r 2 — ri/2 in place of r 2 , while retaining ri as the distance 
between the positron and electron. The Hamiltonian for these asymmetric coordinates is 
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H = - 2 


3, 

4 


K + "vL + - + 


r i 


|R 2 + 5ri| |R 2 


I'll 


(142) 


The final state is given by the p-wave part of the symmetrized product of a plane wave in relative 
coordinates: 


= ^(ri)e lk - R * + 0(r 2 )e‘ kRl ], (143) 

where E — | k 2 . The photodetachment cross section in the velocity form can be written as [75] 

<’ v = ^r\<*f\Qv\9i>\ 2 , (144) 

where the dipole transition operator in the velocity form is 

Qv = 2k • (V P1 + V P2 - V x j, (145) 

and uj is energy of the incident wave. The finite mass of the positron gives a factor of | when 
compared with the photodetachment expression for an infinitely massive atomic ion. The cross 
section in the length form can be written as [75] 


°L = 


2ku>aal 


\<9 f \Q L \9i> 


where the dipole operator in the length form is 


(146) 


Ql = k* (pi + p 2 -x). (147) 

These transition operators can be written in terms of the unsymmetric coordinates: 

Qv — k • (g Vr + 4V ri + 2Vr 2 ). (148) 

and 

Ql — k • (R + ri + ^-^ 2 )* (149) 

Now we represent the initial loosely bound state function in the following form [76,77]: 

g— 7 Rj 

Vi -C— — <f>(rk), for Rj >> r k , (150) 

Hj 

where 7=0.12651775 from our best variational value [68] of the Ps“ binding energy 0.024010113 

Ry. 


C(r) = \/87rre 7r ^(0, r, r). 


(151) 
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Figure 5: Photodetachment cross sections (dashed lines) in the asymptotic approximation for Ps 
and H“ as functions of wavelengths of the incident light. The length and velocity forms of the H 
cross sections are from more elaborate theory of K. L. Bell and A. E. Kingston [78] (solid lines). 


We find C=0. 1856(2) (cf. Ref. [74] for details). We find 

k 3 c 2 

<t v = <tl = (3.8245 x l(T 17 cm 2 ) 2)3 > ( 152 ) 

which can be written as 

<t = (1.32x 10- 18 cm 2 ) (k2 ^ (153) 

We can write the cross section in terms of the wave length of the incident light, 

a = (650 X 10 _18 cm 2 )(^-) 3/2 (l - A-) 3 / 2 , forA < A 0 , (154) 

Ao Ao 

where A 0 =37953.46 A . In Fig. 5 we have plotted the present results [Eq. (154)], compared with 
the corresponding results [77] for H - . In the latter case the asymptotic approximation is seen to 
compare fairly well with the more elaborate calculations of Bell and Kingston [78], and we expect 
the present results to be similarly reliable. 

A scattering calculation for the final state has been carried out by Ward, Humberston and 
McDowell [79] to calculate the photodetachment cross section of Ps - . See the article by Dr. S. J. 
Ward in these Proceedings. 


Muonic Molecules 

There have been speculations of the possibility of realizing useful muonic catalyzed fusion. The 
Born-Oppenheimer approximation has been used traditionally to calculate energy levels of muonic 


31 



Symposium on Atomic & Molecular Physics 


systems. Drachman and I [80], for the first time, using Hylleraas type wave functions for such 
systems, carried out straightforward Rayleigh-Ritz variational calculations and showed that the 
binding energy of even the weekly bound tdfi molecular state of the angular momentum equal to 1 
and vibrational state 1 (J = l,v = 1) can be obtained fairly accuarately. We [81] further studied 
the deexcitation of tdfj, muonic molecule by internal conversion and also calculated the J — 2 
binding energy. We wrote an easy-to-read review article [82] on this subject. I will not discuss 
further our work because Prof. E. A. G. Armour has carried out extensive work on muonic systems 
and has written an article on muonic physics which appears in these Proceedings. 


Polarizabilities of Two-Electron Systems 

We will see in the subsequent section that in order to calculate the Rydberg levels of three- 
electron systems, the interaction between valence electron and core should be known. This inter- 
action can be represented by the potential 


U (*) = -^2 + 


Oil 6/?i — Oi2 


+ 


higher — order terms, 


(155) 


where x is the relative distance between the valence electron and the core, and a\ and c *2 are dipole 
and quadrupole polarizabilities, respectively, and is the first nonadiabatic coefficient. These 
polarizabilities, in the second-order perturbation calculations, are given by the general expression 


Sk,i = £ 

N 


<o\v { \n>< n\v { \o> 
(E n - E 0 f 


(156) 


where |0 > represents the ground state of the core and \N > intermediate states, determined by 
diagonalizing the appropriate Hamiltonian. We can write 


and 


Si,i = a, = £ 


< 0\vi\N >< lV|t7,-|0 > 


N (En - Eo) 


s 2 ,i = /?,• = £ 


N 


< oK|iv >< jv|t7,-|o > 

(E N - E 0 y 


S 3 ,i = 7,- = £ 


N 


< Ol^iliV >< N\vi\0 > 
(E n - Eq) 3 


(157) 

(158) 


(159) 


The potential between the valence electron at a distance x from the nucleus and the core can be 
expanded in the form 


T/ u i , u 2 , 

V — — 2 H 3 + •••• 

X L X J 


(160) 


where 
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and 


Vi = 2 


1 + 


K(Z - 2) 


(w-x), 


(161) 


v 2 = 2 [riP 2 (ri • x) + rlP 2 (r 2 ■ x) - Kw‘ 2 P 2 (w • x)j . (162) 

Here w = ri + r 2 , K — 2/(1 + M), M is the nuclear mass, and Z is the nuclear charge. Drachman 
and I [83] calculated various quantities in Eq. (156) using the pseudostate summation method. We 
used Hylleraas type wave functions for the ground state and intermediate states. The nonlinear 
parameters in intermediate states are optimized by maximizing oq, which according to the varia- 
tional principle, has a lower bound to the exact aq. In Table IX, we give our results for various 
ions. 


Table IX. Adiabatic and nonadiabatic polarizabilities 
for three isoelectronic systems. A(-B)implies AxlO -5 


System 

Li+ 

Be+ 2 

B+ 3 

Oil 

0.192485 

0.052282 

0.019651 

(h 

0.03529 

4.919(-3) 

1.125(-3) 

7i 

6.806 (-3) 

4.847 (-4) 

6.723(-5) 

0^2 

0.11389 

0.01532 

3.427 (-3) 

02 

0.01668 

1.132(-3) 

1.524 (-4) 

72 

2.584 (-3) 

8.819(-5) 

7.136(-6) 


Polarizabilities of He and H 

We now turn to calculations of such quantities for He and H“, where the convergence is rather 
slow for He and even slower for H _ . For He, we [84] used up to 525 terms for the ground state, 364 
terms for the P - wave intermediate states. For D-wave intermediate states, both sd and pp terms 
are required, and therefore we used up to 165 sd terms and up to 56 pp terms in our wave function. 
In Table X, we show our final results for He. 

Table X. Adiabatic and nonadiabatic polarizabilities 



for 4 He. We also give 

the estimated errors. 

k 

i— 1 (Dipole) 

i=2(Quadrupole) 

i 

oi=1.383 241 013 80±10 

02=2.443 372 616 06±20 

2 

/?i =0.707 521 492 749±55 

/J 2 =1.035 440 519 03±20 

3 

71=0.375 538 368 413±21 

72=0.467 345 191 67±23 


Since we include the mass polarization term in the Hamiltonian, it is interesting to see how 
large an effect this has on various polarizabilities. Therefore, we carried out our calculations, using 
the largest expansion lengths, for K = 0 and for finite /i, both for 4 He and 3 He. We present the 
results in Table XI. 
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Table XI. Effect of the finite nuclear mass on the polarizabilities 
of both isotopes of helium. A is the difference between the results 
for finite and infinite mass. 


Quantity 

A( 3 He) 

A( 4 He) 

Oil 

6.4798(-5) 

4.88345(-5) 

Hi 

1.5044 (-5) 

1.1 341 (-5) 

7i 

4.5918(-6) 

3.4638(-6) 

&2 

-2.270054(-3) 

-1.7107096(-3) 

@2 

-8,70263(-4) 

-6.55793(-6) 

72 

-3.64039(-4) 

-2.74316(-4) 


We carried out a similar calculation for H“ for finite K , except that in this case we used 615 
terms in the ground state wave function. The quantities are much larger than in the case of He 
and the results are given in Table XII. 


Table XII. Adiabatic and nonadiabatic polarizabilities 
for H~. We also give the estimated errors. 


k 

i=l (Dipole) 

i=2(Quadrupole) 

1 

«i=206.148 7618±37 

a 2 =7766.79 374±48 

2 

/?i— 1886.699 325±34 

/? 2 =70 155.536 09±25 

3 

7! =20 046.0671±30 

72=711 88.8802±96 


To ckeck the accuracy of our results, we consider two special cases, related to £0,1 and Si^i. 
The first one tests the completeness of pseudostates that we are using without considering the 
ground-state accuracy: 


y: < 0|w ■ x | N >< N |w • x|0 >=< 0|(w • x) 2 |0 > (163) 

N 

Comparison between the left and right sides of Eq. (163) gives a measure of completeness of the 
set N. We give results in Table XIII, along with the difference of the ratio from unity. 


Table XIII. Comparison between left and right side of Eq. (163). 
A— difference of the ratio from unity. 


System 

left hand side 

right hand side 

A 

He 

0.752 552 891 661 

0.752 552 891 770 

1.0(-9) 

H“ 

7.488 423 814 28 

7.488 424 910 34 

1.5(-7) 


The second test is the Thomas-Reiche-Kuhn sum rule, which for the finite-nuclear mass takes the 
form 


Y, < 0|w • x | N >< N |w • x|0 > (Em - E 0 ) — 2 + I< (164) 

N 

Now the extent of the agreement between the left and right sides of Eq. (164) measures both the 
accuracy of the ground state |0 > and the completeness of intermediate states | N >. We find 
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the ratio of the left side to the right side is 0.999 999 87 for He, and the corresponding ratio for 
H _ is 0.999 9787, respectively. All these tests show that H“ is a more sensitive system than He. 
Nevertheless, our results are accurate to the accuracy given. 


Polarizabilities of and D^ 

Having calculated polarizabilities of two-electron sytems, it must appear that there should be 
no difficulty in calculating polarizabilities of molecular ions H^ and D^, the only difference being 
that two light particles have been replaced by two heavy particles and the mass-polarization term 
becomes as important as other kinetic energy terms in the Hamiltonian. Since, with a hammer 
in hand the whole world looks like a nail, I just went ahead with the calculations for H^ and 
D^" systems using Hylleraas type wave functions, with the appropriate change of masses in the 
Hamiltonian. But I found no expansion length in these functions was adequate enough for ap- 
proaching the known results obtained using the Born-Oppenheimer approximation, leaving aside 
to surpass them. The problem remained dormant till Drachman heard from Steve Lundeen about 
his experimental results on high Rydberg states of H 2 . It is possible to extract accurate properties 
like quadrupole moment and the parallel and perpendicular polarizabilities of the molecular ion 
Hj from these results [85,86]. The Born-Oppenheimer results [87] disagreed with the experimental 
results. We tried to resurrect our old work and tried various extrapolations, but had no success. 
Drachman came up with an excellent idea: simply raise r \2 (the interparticle distance between 
similar particles) to a power close to 30 or so and choose the nonlinear parameter a in e~ ari2 equal 
to half of that power of r 12 in the Hylleraas functions being used. This was a miracle in the sense 
that nearly six terms in the Hylleraas expansion now were equivalent to hundreds of terms in the 
earlier expansion! It is easy to understand, after the fact, that now two protons stayed clamped 
at their respective positions whereas they enjoyed the same freedom as electrons with the usual 
generalized Hylleraas functions. This was just what was needed to get excellent results without 
making Born-Oppenheimer like approximations. The expression for polarizability [88] is given by 


_ 4 1 < 0|£A + %b\N > | 2 3 

ai M 3 (l + v ) 2 'V En - Eq a ° 


(165) 


where za = ta • £ and zb = r# • €, € being the direction of the external electric field which can be 
considered to be in the z direction, ta and r# are the distances of the electron from protons A and 
B, respectively. The reduced mass fi — Mj (M + 1). It should be pointed out that we are treating all 
the particles on an equal footing and we do not refer to any special ’’molecular” quantum numbers. 
Thus we are not interested here in the ’’axial” or ’’transverse” polarizabilities that appear in the 
Born-Oppenheimer approximation. 

With these modified Hylleraas functions, we obtained ground state energy £'(H2")=-1.194 277 
909 Ry, differing by only about 2.2xl0 -7 Ry from the accurate value [89]. We calculate the 
intermediate P states also using high powers of ri 2 (cf. [88] for details). Our final results for H 2 " 
and D 2 " are given in Table XIV. 
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Table XIV. Polarizabilities of and obtained in different ways, 
both theoretically and experimentally. Quantities in perentheses 
are errors in the last decimal place given. 


Method 

ai(H+) 

«i(Dj) 

Experiment 

3.1681(7) 

3.0712(7) 

Born-Oppenheimer b 

3.1713 

3.0731 

Modified generalized 

3.1680 

3.0671 

Hylleraas functions 
Finite element 0 

3.1682(4) 

3.0714(4) 


a Jacobson et al. [86]. ^Bishop and Lam [87]. c Shertzer and Greene [90]. 

It is obvious that the results for are not as close to the experimental results as in the case 
of H^, showing that even the modified wave functions have limitations: as the nuclei get heavier 
it becomes more difficult to generate a well-enough localized two-nucleus part of the trial wave 
function. Clearly, the method fails as the nuclei become infinity heavy. 

Dr. J. Shertzer discusses in these Proceedings the results obtained by the finite-element method. 


Polarizabilities of HD+ 

After we had studied the polarizability of the homonuclear molecular ions it should have been 
simple to extend our methods to the heteronuclear ion HD+. Janine Shertzer reminded Drachman 
that in this case there would be dipole coupling between rotational levels with J— 0 and J— 1, which 
was not possible in the earlier cases because of symmetry. In addition she pointed out that the 
lowest-lying J— 1 level is so close to the ground state that one would expect a very high value 
of polarizability, since the denominator of the second-order perturbation sum would be so small. 
Drachman and I [91] carried out calculations to see if this prediction was borne out. 

The dipole polarizability oq is given by the second-order perturbation theory: 


Oil 


< 0\d\N >< AT|rf|0 > 3 
2. s Z7__ IT « 0 - 


N 


En - E 0 


(166) 


The dipole operator d depends on masses of the nuclei and is given by 

d = 2 


M d + 1 Mp + 1 

r i + —n — r 2 


Mj 


Mj 


(167) 


where the unit vector e is in the direction of the applied electric field, Mp and Mp are the masses 
of the proton and deuteron nuclei and 


Mp = Mp + Mp + 1. (168) 

The calculation of the dipole polarizability of HD + is similar to that of and D^ : we treat all 
three particles on an equal footing and do not refer to any special ’’molecular” quantum numbers. 
We use the modified generalized Hylleraas type wave functions, i.e., we use very high powers of n in 
r i 2 € ~ ari2 an d ot is of the order of n/ 2. The energy eigenvalues using the Rayleigh-Ritz variational 
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principle are J5o=- 1.195 795 889, using 560 terms, and E p =- 1.195 372 602 Ry, using 728 terms, for 
the ground S state and the lowest P state, respectively. These compare well with i?o=-1.195 795 
931 and E p —~ 1.195 396 2560 Ry obtained by Moss [92]. The energies of these two states are so 
close to each other that they are almost degenerate states of opposite parity. This is the important 
difference between the homonuclei Hj and , and heteronuclei HD+ molecules. 

The polarizability due to the lowest P term alone in Eq. (166) is 392.0814 «q and the sum over 
the remaining intermediate states contributes a s =3.2076aQ. The total polarizability, aq=395.289aQ, 
is very much larger than might be expected from the results [88] in and D^. This can be traced 
to the fact, as indicated above, that there is a coupling between the two lowest S and P states 
due to their closeness in energy. If we exclude the ground state, we can calculate the polarizability 
of the lowest P state by summing over all the intermediate S states. We find this result to be 
a p =2.03008ao. 

But second-order perturbation theory can only be legitimately carried out if the perturbation 
is small compared to the spacing of the unperturbed energies, and this condition is hard to satisfy 
in the present case. If the perturbation is due to a unit charge it must lie further from the ion 
than about 32«o for the perturbation theory to be valid. What happens for larger fields or smaller 
distances? Clearly, for these larger fields the opposite parity levels act essentially as degenerate 
levels, so the techniques of degenerate perturbation theory must be applied. This leads to an energy 
shift that is linear in the electric field, falling like R~ 2 rather than R~ 4 as we expect for ordinary 
polarizability. This is usually described as the effect of a permanent dipole moment. We then went 
on to diagonalize the perturbation matrix in various approximations and were able to give a good 
description of the energy shift as a function of R. 

In observing the Rydberg states of HD, the effective polarizability should be a*=3.2076ao, which 
is close to the polarizabilities of and Dj [88]. With this value, we have calculated energy shifts 

AjVL = -<*s J (169) 

using hydrogenic wave functions. We give these energy shifts A ml of a series of Rydberg levels 
due to the polarizability of the IID + ionic core omitting the lowest rotationally excited level in the 
summation over intermediate states. 


Table XV. Energy shifts A nl 
of a series of Rydberg levels. 


N 

L 

-Ajvl(MHz) 

9 

4 

11501 


5 

3943.6 


6 

1594.9 


7 

720.08 


8 

350.41 

10 

4 

8527.4 


5 

2951.8 


6 

1208.8 


7 

554.8 


8 

275.89 


9 

145.2 
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Figure 6: Diagrams symbolizing the nonrelativistic calculation of a(u>), with energies of each state 
indicated. 


Optical Properties of He Including Relativistic Corrections 

An interesting application of polarizabilities is the calculation, which Drachman and I [93] 
carried out, of the index of refraction n of He from which the Verdet constant F, which measures 
the rotation of the plane of polarization in the Faraday effect, can be obtained. The rotation 0, 
measured in degrees, is given by 

0 = VBL , (170) 

where B is the magnetic field and L is the length traversed by light. If the static field is replaced 
by an oscillating field (an electromagnetic wave), it is possible to define a frequency-dependent 
polarizability a (u>) from which the index of refraction can be obtained (cf. Fig 6). The expression 
for the polarizability is given by 


a 


M - ~Y^<o\z\n><o\z\n> 


1 


+ 


N L Ej\f — ( Eo + u>) (En + 2u>) — ( Eo + u>) J 

< 0|Z|iV >< N\Z\0 > (E n - Eo) 


^ (En ~ Eo ) 2 - oj 2 

where N = p and the dipole operator Z is given by 

Z = 2(z\ + Z 2 ). 

We can define a set of ’’generalized dipole polarizabilities” as follows 


(171) 


[«!)/ l>7l> X (E n - Eo)( 1-2, 3, 4, 5,6,7)’ 


(172) 


(173) 
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As indicated earlier, we use Hylleraas basis sets for the ground state as well as for pseudostates N 
and the results for various quantities are given in Table XVI. 


Table XVI. Various quantities 
for 4 He. 


Oil 

1.383 241 01 

(h 

0.707 521 493 

7i 

0.385 538 368 

Sx. 

0.218 735 026 


0.127 538 649 

& 

0.075 827 657 

m 

0.045 731 135 


The denominator in Eq. (171) can be expanded in powers of u/(Ejy — Eo) <0.2 for wavelengths of 
visible light: 


a{uj) = g?i + 71 (jJ 2 + ciu ; 4 + rjxuj 6 + ... (174) 

The relativistic corrections are calculated by using the Breit-Pauli relativistic Hamiltonian which 
has the following form: 

B = a 2 + V 4 ] + 27r[«(n) + S( r 2 )] + 2^(r 12 ) + ~[V l • V 2 + r 12 • (r i2 • Vi)V 2 ]} . (175) 

This requires a third-order-perturbation theory, with B retained to first order and Z to the second 
order. Up to u> 4 the expansion corresponding to Eq. (174) gives the following numerical result 

Aa(«) = -8.0029 X 10~ 5 - 8.1516 X 10~ 6 u> 2 + 3.006 X 10“ 7 u; 4 , (176) 

and the relativistic expression for the frquency-dependend polarizability of helium becomes 

ctreiM = 1.383160981 + 0.385530216u; 2 + 0.12753895w 4 + 0.045731 14u> 6 . (177) 

The last term Eq. (177) has not been modified from its nonrelativistic value, since the effect of 
relativity here would be absolutely negligible. In order to write the above expression in terms 
wavelength, we use 


47ra 0 911.39198 

a{ 1 - Kf 2) A ~ A 


(178) 


where /i'=2.741 493x10 4 for 4 IIe and the wavelength is in A units. We obtain the expression for 
polarizability in terms of wavelength in the form: 


a re i (A) = 1.383729930 + 


3.204546 X 10 5 8.822907 X 10 10 2.624092 x 10 16 


A 2 


+ 


A 4 


+ 


A 6 


(179) 
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Optical measurements usually give n — 1, where n is close to 1 for helium. The higher-order 
corrections in the relation between the polarizability (a single atom property) and the index of 
refraction, due to the effects in the medium, are accounted for by the Lorenz-Lorentz equation: 

i 4-7T 

— — = —N L a 3 0 = 16.67718 X l(T 6 a r ei(A) = 2 (180) 

Tl 1 A O 

where Nl is Loschmid’s number. Solving for n to second order in z, we find 


3 3 y2 

n- 1 ~ -z + -- = 34.61527 x 10" 6 
2 8 


8.016511 2.207154 x 10 6 6.564503 x 10 11 , 

+ 77 — + 77 + 77? . (181) 


A 2 


A 4 


A 6 


The Verdet constant is given by 


e dn _ 1.616813 X 10 7 
2 me 2 d\ A 2 


1 + 


5.506521 x 10 5 
A 2 


+ 


2.456618 x 10 11 
A 4 


/i min/oerem, (182) 


where min refers to rotation. We get from the above equation V=0.661 22±0. 00025 at A=5000A 
versus the experimental value 0.637 of Leonard [94], while at A=8000A we get V=0.254 816±0.000 
015 compared with the experimental value 0.246 (all in units of n min/oer cm). This shows 
that there is a significant discrepancy here that is larger than the relativistic effects that we have 
considered. The experimental results do not appear to be very accurate and new measurements are 
required of the refractive index n of helium to have a better comparison between our theoretical 
results and measurements. 

Another Way to Calculate Lamb Shift 

One of the most difficult parts in two-electron Lamb shift calculations is the Bethe logarithm 
given by 


ln(K) = 


< 0|Vlrc >< n\V\0 > (En - £o) 3 ln(E n - E 0 ) 
J2n < 0|V|n >< n\V\0 > {E n - E 0 ) 3 
N_ 

D ’ 


(183) 


where | n > are L=1 eigenstates, both bound and continuum, of the Hamiltonian describing the 
two-electron system. Instead, we use the pseudostate summation method. The interaction V is 
given by 


V = zi + z 2 = ricos#i + r 2 cos0 2 - (184) 

Having used pseudostates in the polarizability calculations, it should have been easy to calculate 
ln(A) given by the above expression in the length form. But the result obtained using either the 
length, or velocity, or acceleration form never seemed to approach the known results, no matter how 
large the psoudostate expansion was in the above expression. The basic reasons are the power of 
(E n — E 0 ) which is 3 instead of -1, as in the polarizability expression, and the presence of \n{E n — Eo) 
in the numerator, which makes matters worse. Combined together, the convergence of the results 
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becomes extremely slow. Out of frustration, I sought Drachman’s help. To solve this problem, he 
came up with several tricks which I prefer to call good insight. We could then obtain results as 
accurate as available in the literature. I have used a slightly different notation in this section for 
the pseudostates in order to conform with our published paper [95]. 

The Denominator 

We transform the denomenator D in Eq. (183) in a form which does not have an intermediate 
sum | n > by using the commutation relation 

< n\[H, V]\0 >=< n\V\0 > ( E n - E 0 ) (185) 

three times, in order to remove (E n — Eq): 


D = |(2 + A’) 2 {< 0|[ff, (V x + V 2 )] • (Vx + V 2 )|0 > - < 0|(Vx + V 2 ) • [H, (V x + V 2 )]|0 >} . 

(186) 

Since the potential term in the Hamiltonian does not commute with Vi + V2, the above expression 
can be put in the form 


<0|5( ri ) + 5(r 2 )|0>, (187) 

where closure over the intermediate states | n > has been invoked and Poisson’s equation has been 
used to introduce S functions. As before, K — 2/(M + 1), where M is the mass of the nucleus 
in units of the electron mass. The terms of the order K 2 have been dropped; within this error 
the expression is exact. It does depend on the accuracy of the ground state wave function which 
consists here of 525 terms. Therefore, the error due to the accuracy of the ground state wave 
function is negligible. 


Table XVII. ln(A' ) as a function of 7 and Pekeris number £l p . 


ftp 

7 

D a cc 

N acc 

ln(A) 

3 

3.644 529 

94.114 569 

301.053 261 

3.198 796 

4 

4.080 394 

101.216 416 

346.195 202 

3.420 346 

5 

4.425 096 

104.957 668 

373.513 205 

3.558 703 

6 

4.849 542 

108.413 460 

399.676 322 

3.686 593 

7 

5:179 484 

110.403 345 

416.057 508 

3.768 523 

8 

5.588 226 

112.320 411 

432.319 315 

3.848 983 

9 

5.888 226 

113.508 976 

442.840 189 

3.901 367 

10 

6.305 080 

114.676 534 

453.744 402 

3.956 373 

11 

6.623 156 

115.446 146 

461.125 678 

3.994 293 

12 

7.004 813 

116.207 693 

468.586 094 

4.032 316 

13 

7.3 

116.735 692 

473.876 067 

4.059393 


Because of the In (E n — E 0 ) term, it is not possible to reduce N to a similar form. However, 
we expect the convergence of pseudostates will be similar for N and D. We use the commutation 
relation [Eq. (185)] twice for each matrix element of Eq. (183) to obtain the acceleration form 
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Figure 7: Behavior of D acc for helium as a function of the nonlinear parameter 7 . The crosses are 
for and the plus signs are for fi=4. The dotted line is the ’’exact” value, obtained from Eq. 

(187). 


D 


acc 


-8Z 2 (1 + K)J2 

n 


< 0\U\n >< n\U\0 > 
Eq ~ E n 


(188) 


where 


U = 


COS01 COS02 
“ 72 “ + 


1 


t*2 


(189) 


This expression for D acc looks like the second-order energy shift induced by the potential U , and 
it has a variational bound [95]. We can, therefore, choose the nonlinear parameter 7 in the pseu- 
dostates to maximize D acc in Eq. (188) (cf. Fig. 7). 

As the number of \n > is increased the D acc approaches D ex which for He from Eq. (187) is 
121.335 143. We calculate N for the same 7 and the same number of | n > as for D acc . We give 
values of ln(/i ) as a function of the Pekeris number in Table XVII. 

Extending The Upper Limit 

We assume that the contribution of each term in the pseudostate summation is exact and the 
remaining error in the total is due to the fact that that the sum does not extend to infinity. We 
use the method of Dalgarno and Stewart [96] to account for the remainder of the sums beyond the 
highest pseudoenergy. For high energies they used the following simplified form: 


1 ZtA 1 / 2 

®n,k(ri,r 2 ) = — [u„(r 1 )e ikr2 + tin(r 2 )e' k ‘ ri ], (190) 

to represent the singly ionized states in the expression that replaces the discrete one in Eq. (188): 
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F n (e) 


16Z 2 (1 + K) 
3 (I n + *) 


< ®o| 



|®n,k > 


(191) 


Here e = k 2 , I n is the ionization-excitation potential of two-electron systems, and u n is the wave 
function of the one-electron system left after single ionization. Without loss of accuracy, we fitted 
F n (e) to the form [A n + (B n /A;)tan -1 (Cn/k)]€~ 3 / 2 and included s states up to n= 4. Higher states 
were included approximately by an exression falling like 1/n 4 . We adjust €q so that 


roo 

D ex — D acc -\- / de F n (e). (192) 

n 

The critical step is now to correct N using the same value of e 0 as for the denominator, modifying 
the integral by the inclusion of the appropriate logarithmic factor. In Table XVIII, to obtain the 
exact value of the denominator we give the required value of €o, the corrected value of the numerator 
and the improved value of In (A'); the convergence with Q p is significantly improved. 


Table XVIII. The €q, corrected values of iV corr , and In (A"). 


D, p 

£0 

N corr 

\n(I<) 

3 

361.0613 

520.267 518 

4.287 855 

4 

724.0096 

521.248 382 

4.295 939 

5 

1143.441 

523.110 217 

4.311 284 

6 

1912.977 

524.085 779 

4.319 324 

7 

2734.433 

525.086 045 

4.327 568 

8 

4108.808 

525.799 439 

4.333 447 

9 

5524.324 

526.334 777 

4.337 859 

10 

7729.503 

526.911 785 

4.342 614 

11 

9964.945 

527.307 234 

4.345 874 

12 

13 253.52 

527.650 177 

4.348 700 

13 

16 564.86 

527.870 611 

4.350 517 


Extrapolation and Results 

We extrapolate our results by using the deviation AD of the denominator from its exact value: 


AD = D ex - D acc (193) 

This quantity approaches zero as D p increases. We plot the uncorrected and corrected In (A') (cf. 
Fig. 8). 

The slope of the line fitted to the corrected points is 12 times smaller than the uncorrected 
slope, indicating the improved convergence we expected. The extrapolated result for helium is 
ln(A)=4.367 58(46). 

In Table XIX we show the values of In (A'), obtained by this method, for a range of atomic 
numbers Z and compare them with those obtained from the approximate expression of Goldman 
and Drake [97]. Their approximate expression has been obtained using an expansion in 1/Z: 

ln(K) = ln[19.7692669(Z - 0.00615) 2 ]. (194) 
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Figure 8: ln(/i ) for helium as a function of AD. The plus signs show the uncorrected values, using 
the data from Table XVII, and the crosses are the corrected values from Table XVIII. Linear fits 
have been made in both cases, and the improved convergence of the corrected values is clear. 


Table XIX. Results for a series of two-electron systems. 


System 

Z 

In (A') 

ln(A') a 

He 

2 

4.367 578 

4.364 263 

Li+ 

3 

5.177 763 

5.177 249 

Be+ 2 

4 

5.753 615 

5.753 640 

Ne+ 8 

10 

7.586 072 

7.588 068 


a From Eq. (14) of Ref. [97] which is the same as Eq. (194) given above. 


A good fit to our results is the form: 


ln(K) = ln[19.705541(Z+ 1.35 X 1(T 5 ) 2 ]. (195) 

The present results may represent the two-electron ground state Bethe logarithm well over the 
range 2< Z <10. 


ln(R~) for Atomic Hydrogen 

It is an obviously interesting question to see how well the above method works for the simpler 
and better known one-electron (atomic hydrogen) case. Here the ground state wave functtion is 
known exactly and the pseudostates are of a simple one-electron form: 




7T 1 / 2 ’ 


V n = e-T r P l (0)Y t C j r n . 

3 = 1 


(196) 
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Figure 9: ln(K) for hydrogen as a function of AD. The alternating convergence pattern discussed 
in the text leads to two distinct lines: crosses are for even values and plus signs are for odd values 
of 


The denominator can be evaluated exactly: D=16/3. All other steps described above can be 
applied to the one-electron system. 

In Fig. 9 we show the unexpected results. In place of the linear relation found for two-electron 
systems, we obtain an alternating convergence pattern. As increases from an odd value to the 
next even value the numerator increases significantly while the denominator remains unchanged 
to 7 or 8 significant figures. With the next increase in the increase in D resumes. This effect 
produces two distict lines. We were able to carry the calculations up to fi p =22, at which point 
AD=0.057. The extrapolations to AD=0 give ln(K)=2.987 125 (from odd values of fi p ) and 
ln(A”)=2.978 329 (from even values). Combining the two results we can report a ’’best” value of 
ln(A r )=2.9827 ±0.0044 which should be compared to the accurate value of Haywood and Morgan 
III [98] ln(A')=2.984 129; our generous error does include this value. 

We have not been able to understand this irregular convergence and this remains an interesting 
unanswered question. 


Rydberg States of Li 

Traditionally, eigenvalues and eigenfunctions are calculated by the use of the Rayleigh-Ritz 
variational principle. This procedure has been carried out for states with high quantum numbers, 
N and L, as well [99]. The disadvantage is eigenvalues of all the states below the state of interest 
have to be calculated. Drachman realized that when N is large, say 10, the outer electron is so far 
away [100 ciq compared to ao for the core electrons] from the spherically symmetric core that it does 
not have much electron-electron correlations, the type taken into account by the Hylleraas functions, 
with the core electrons. Even the exchange is not important, and most of the of correlations can 
be taken into account by considering only the long-range interactions. On this basis, he developed 
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a formalism [100] to calculate energies of states of interest only, using the Feshbach projection 
operator technique. Since the exchange can be ignored Eq. (71) simplifies to 

P = Pi = \<f>o><M- ( 197 ) 

The formalism is rather complicated and I give the final result for the effective potential seen by 
the outer electron at a distance x from the nucleus: 


TT , x «i , 6g/?i - « 2 , $ + 169 2 7i/5 

U \ x ) = -3T + ZZ + It 


+ 


-a 3 + 15 qf3 2 - e + a v 3i - 72^ 2 7i[1 + L(L + 1)/10] 


(198) 


I alluded to this form in Eq. (155), where I indicated the importance of polarizabilities in the 
Rydberg states. The core coefficients are described below: 


c . = ~2i+2 0|v;-|» >< n|V|0 > 

k ’ i - X ^ (E n -Eo) k 


(199) 


Here, as mentioned before, a* = and 7 ,• = S^-. The third-order polarizability has 

the form: 


6 = 


* 7 £ 

n,m 


< 0|V|n >< n\Vj\m >< ra|Vt|0 > 
(■ E n - E 0 )(E m - E 0 ) 


(200) 


The values that ( ijk ) can take are all the permutations of [112]. Finally, the fourth-order hyper- 
polarizabilty involving the dipole terms has the form: 


_ 8 r' < 0|Fi|n >< n|Vj|m >< m|Vi|p >< p\V\\0 > 

£ (E n -Eo)(E m -Eo)(E p -E 0 ) ' 1 j 

Drachman [100] developed the formalism for helium atoms and showed that he could obtain 
the same result as Drake [99] did for the eigenvalues, e.g., iV=10 and £=6 his result £=-105.829 
80±0.00014 compares very well with -105.829 683 489 MHz of Drake (cf. Table III of Ref. [100]). 

Up to now I have not given the interaction potentials Vi, V 2 , and V 3 . We are now interested in 
the Rydberg states of Li where the spherical core is He-like and therefore, I give below expressions 
[101] for three-electron systems: 


V, 


2q 2 w 


2+ (Z- l)K 


V 2 = 


2q 3 


2+ I< 

r'lPii ri -x) + r 2 P 2 (r 2 • x) 


2 u tt . _ ^ , „2 __ 4K ± ZI< 2 -.2 

(2 + K) 2 


Pi(w-x), 

A-w 2 F 2 ( w -x) 


Vs 


2q 4 [ 2 — £' 3 3 .... 4K-2K 2 -ZK 3 3n ,„ 

-[r 1 P 3 (ri • x) + rfP 3 (r 2 • x)] w 3 P 3 {w • x) 


2 + K 


(2 + A') 3 


(202) 

(203) 

(204) 
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where w = ri + r2 and q ~1+K 2 /2. In Table XX we give the various quantities for 7 Li+ and 
6 Li + , which have been calculated by using the appropriate Hylleraas type wavefunctions for the 
ground state and intermediate pseudostates of angular momentum L— 1, 2, and 3. Again we use 
the method of pseudostate summation. 


Table XX. Best values of the core parameters used in 
constructing the effective potential. The upper 
entry is for 6 Li and the lower one is for 7 Li. 


Oi 1 


7i 

Oi2 

0.192 490 771 
0.192 485 410 

0.035 286 879 
0.035 286 017 

0.006 806 377 
0.006 806 227 

0.113 825 934 
0.113 834 685 

02 

72 

5 

€ 

0.016 670 328 
0.016 671 511 

0.168 351 237 
0.168 362 339 

0.121 .337 559 
0.121 345 411 

0.027 039 600 
0.027 038 073 


From the core parameters of Table XX the effective potential of Eq. (198) can be constructed 
explicitly. Numerically, this potential is 

TTf 0.19248540 0.097881 0.143125 0.428584 + 0.049005L(L + 1) , onK ^ 

U[X >- ^ + X 6 +— p W 

for the case of 7 Li. We use the following expression for the energy shift A (NL) away from the 
unperturbed energy - R/N 2 


A (NL) = R 


(U 4 + Ue) + \[U 7 + U s ]) ± \[U 7 + Us]} , 


(206) 


where the reduced Rydberg R= 3.289 584 678 xlO 9 MHz for 7 Li, and where Uk is the expectation 
value of that term in U(x) which goes like x ~ k . Since the unperturbed wavefunction of the outer 
electron is purely hydrogenic, these expectation values can be evaluated analytically and exactly. 
Some results are given in Table XXI for N — 10 and various values of L (cf. Table II of Ref. [101] 
for values for other N and L). 


Table XXI. Level shifts (in MHz) for 7 Li for iV=10 
due to the effective potential terms Uk- The total 
and error are obtained as described in Eq. (206). 


L 

U . 4 

U 6 

U 7 + U 8 

Total 

Error 

4 

-511.674 

2.332 

-0.229 

-509.457 

0.115 

5 

-177.1181 

0.3007 

-0.0022 

-176.8185 

0.0011 

6 

-72.5341 

0.0541 

0.0005 

-72.4798 

0.0003 

7 

-33.29026 

0.01190 

0.00014 

-33.27829 

0.00007 

8 

-16.554171 

0.002929 

0.000033 

-16.51226 

0.000017 

9 

-8.712722 

0.000738 

0.000007 

-8.711981 

0.000004 


Since Us has been included, a.\fx 4 should also be included to the second order. That is 
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/ 


a 2 = E 

N 


< NLlotx/x^N'L >< N'L\a 1 /x 4 \NL > 
En - En> 


( 207 ) 


Here the intermediate states N'L are the hydrogenic states of the outer electron. This quantity 
has been calculated for He [100], A 2 for 7 Li can be obtained by scaling the results for He; A 2 being 
proportional to a 2 . The results, given in Table XXII, are again shown only for the N=10 and 
various L. The large relativistic correction due to increase in mass with velocity has been discussed 
by Bethe and Salpeter [102] and some details are given in [101], as well. We obtained the following 
expression for the leading relativistic correction: 


A re l — 


a 2 R 

IP 


JLf! -—] — 

4 iV L 6 J 2L + lj ’ 


which is accurate up to order K 2 . Again some results are given in Table XXII. 


(208) 


Table XXII. Second-order corrections, leading realtivistic 
corrections, and the final total shifts for 7 Li in MHz for N =10. 


L 

a 2 

A re i 

Total-)- A 2 + Arel 

4 

-0.097 

-25.790 

-535.343i0.115 

5 

-0.0090 

-18.7122 

-195.5397±0.0011 

6 

-0.0012 

-13.8122 

-86.2932i0.0003 

7 

-0.0022 

-10.218 88 

-43.497 39i0.000 07 

8 

-0.000 046 

-7.471 035 

-24.022 31i0.000 02 

9 

-0.000 Oil 

-5.301 687 

-14.013 679i0.000 0004 


The fine-structure splitting for the N = 10 manifold are shown in Table XXIII, where they are 
compared with the accurate measurements of Rothery et al. [103]. 


Table XXIII. Fine-structure splitting for the N — 10 manifold of 7 Li . 


Interval 

Enegy shift (MHz) 

Experiment [103] (MHz) 

10G-10if 

339.80i0.011 

339.7186i0.0031 

10/7-10/ 

109.246 6i0.001 1 

109.2140i0.0047 

10/- 10/C 

42.795 8i0.000 3 


10A-10L 

19.475 08i0.000 07 


10L-10M 

10.008 63i0.000 02 



Considering the accuracy of the measurements, the agreement between theory and experiment 
is not very good. Drachman and I [104] extended our work [101] to included corrections to the 
third order to improve the agreement. 

Relativistic Correction to the Polarization Potential 

We wish to compute a correction to the energy shift of a Rydberg level due to the Breit-Pauli 
relativistic Hamiltonian of order a 2 . Thus we must carry out a third-order perturbation calculation 
which has the form 


T ,/_\ ^<0\h\n><n\h\m><m\h\0> ^ < 0\h\n >< n\h\0 > 

ViM = L (£ 0 2- < 0|*l» >■ < 20SI ) 


(Eo - E n y 
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where h^II^ip + 7/bp • The expression for the Breit-Pauli Hamiltonian Hbp = B has been given in 
Eq. (175) and 


Hd\ P - -o(ri + r 2 ) • x. 


Keeping terms to order a 2 results in the following: 


( 210 ) 


V s ° 


= E 


< 0|i?dip|ft >< n\HBp\m >< m|Hdip|0 > o < 0|#di P |ft >< n\Hdi p \m >< to|.J7bp|0 >" 

I ^ 


(Eo - E u )(Eq - E m ) 

E < >< »i^- liP )° > < 0|ffBp|0 > 

n \Eo E n ) 

Acq 


(Eq - E n )(Eo — E m ) 


(211) 


where we have taken account of the facts that the ground state has angular momentum L— 0, the 
excited states n, m are L — 0 or 1, and i/g P is rotationally invariant. Our best-converged result, 
obtained with 161 terms with L— 0 and 165 terms with L= 1 is Aai =-4.518x 10” 5 . 

Retardation Corrections 

As indicated earlier, the Rydberg electron is at a distance from the nucleus much greater 
compared to the radius of the core. When the distance is greater than 137ao, the interaction is no 
longer purely Coulomb in character. This is because the delay due to the finite light propagation 
time between the core and the outer electron is comparable to the characteristic time t = ao/v. 
This retardation (or Casimir effect) brings in a new type of term [105] in the effective potential 
acting on the Rydberg electron that falls off like x“ 5 . Au, Feinberg, and Sucher [106] have given 
the following expressions for the modification of the effective potential producing the energy shift 
in the state (1 s 2 NL): 


/\ NL 

A Ret 


16 


E 


|< n\zi + z 2 l0 >| 

(E n ~ Eo) 2 


( 212 ) 


In 



2znt [3 - 5 z 2 + z* 4 - (6 z n - 2 Zn)t] 


37T 

~2 


(213) 


where z n —^ax(E n — Eo)- The evaluation of this correction is rather complicated. Nevertheless, 
we have evaluated the retardation corrections from N=b to IV— 21. We give results for the iV=10 
manifold only in Table XXIV. 


Table XXIV. Retardation corrections for lithium in MHz for iV=10 manifold. 


L—9 

L—8 

L—l 

i II 
! ^ 

T= 5 

L=4 

0.000646286 

0.00142682 

0.00325464 

0.00790948 

0.0212477 

0.0653658 
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Lamb-Shift Corrections 

Since we are interested in the T-dependent fine-structure splitting of the Rydberg levels only the 
change in the Lamb shift of the two-electron core due its interaction with the outer electron needs to 
be calculated. The main parts of the Lamb shift (mass renormalization, vacuum polarization, and 
radiative corrections to the magnetic moment) can be written in terms of (< S(ri) >+< 5(r2) >), 
and it is necessary to calculate the dependence of these 8 functions on the state of the outer electron. 
Following Goldman and Drake [107] we can write the expression for the two-electron Lamb shift as 


A NL _ 8 7 2 
A Lamb — 


19 

-21na + — - InK 
30 


(< <H r l) > + < <H r 2) >)• 


(214) 


Here the 8 functions refer to the two core electrons but are influenced by the outer electron. This 
correction is proportional to the expectation value of 1/x 4 and behaves like another correction to 
the dipole polarizability of the two-electron core. We can write the relativistic corrections as 


A N l = [1-486 x 10 5 - 3.103 x 10 4 ]/lj\ + A^, (215) 

\ x / NL 

where the quantities in the square bracket are the coefficients (in MHz) of the relativistic po- 
larizability and Lamb-shift corrections, respectively. In Table XXV we show the three types of 
corrections for the experimentally interesting N = 10 and their total in MHz. 


Table XXV. Relativistic, polarizablity, retardation, Lamb shift corrections 
and the total uncorrected interval from Table XXII, in MHz, for N = 10 manifold. 


L 

Uncorrected 

shift 

Relativistic 

polarizability 

Retardation 

Lamb shift 

Total 

4 

-535.343±0.115 

0.1201 

0.0654 

-0.0252 

-535.183i0.115 

5 

-195.5397±0.0011 

0.0416 

0.0212 

-0.0087 

-195.486i0.0011 

6 

-86.2932i0.0003 

0.0170 

0.0079 

-0.0036 

-86.2719i0.0003 

7 

-43.49739±0.00007 

0.0078 

0.0033 

-0.0016 

-43.4879i0.0001 

8 

-24. 02231±0. 00002 

0.0039 

0.0014 

-0.0008 

-24.0178 

9 

-14.013679±0.000004 

0.0020 

0.0006 

-0.0004 

-14.0115 


Finally, in Table XXVI we compare the experimental fine-structure intervals for lithium [103] 
with the theoretical totals including the uncorrected values and the three small corrections. It is 
clear that there is better agreement when the small corrections are included. However, higher-order 
corrections and more measurements seem to be warranted. 


Table XXVI. Comparison of level differences for lithium, 
in MHz, between theory and experiment [103]. 


Interval 

Experiment-Theory 

Standard deviation 

lOG-10/f 

0.02 

0.11 

10H-10I 

0.0003 

0.0048 


We [108] have carried out a similar calculation for CIV, O VI, and Ne VIII. But at present there 
are no measurements accurate enough to compare with our calculations. 
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Even though there are other interesting topics on which Dr. Drachman and I have worked 
on together, I stop here. Instead, I mention below two more topics: one with others and one by 
himself. 


Positronium-hydrogen scattering resonances 

Drachman once tried to insert the statement ’’Nobody likes a smooth cross-section” in an article 
on resonances, but the referee vetoed it. Nevertheless, there is much more interest in scattering 
resonances than in bland nonresonant behavior. The Ps-H system is rich in interesting physics: 
There is one particle-stable state called positronium hydride [PsH] with an energy of about 1 eV 
below the free Ps+H threshold. Using both the stabilization method and the complex rotation 
method Drachman and Houston [109] found an s-wave resonance in elastic Ps-H scattering at 
about 4.5 eV. 

Since interesting physics usually emerges from the analysis of resonances, it was of interest to 
understand the mechanism producing this resonance. At first, it was thought that it was due to 
some threshold process like that producing resonances in electron-hydrogen scattering below the 
n—2 threshold, but the position obtained was not close enough to a threshold to make this plausible. 
The best explanation describes it as a Feshbach resonance [110] in which the closed-channel part 
of the scattering function is the re-arranged system e + -H _ , corresponding to perturbed hydrogenic 
bound states. Because of the long-range Coulomb potential between the positron and the negative 
hydrogen ion one can predict that there should be an infinite series of resonances, of which the one 
found by Houston and Drachman [109] is just the first. Since the hydrogen ion exists only in the 
singlet spin state these resonances should not occur in the triplet state. No reliable calculation has 
found triplet resonances. 

There are more subtleties in this system, including some problems with the low-lying resonances 
expected for L >0 [111] and these have been examined very recently by De Rienzi and Drachman 
[112,113]. More about this interesting system can be found in the presentation by Dr. H. R. J. 
Walters in this volume. 


Hyperfine Splitting in Muonic Helium 

Huang and Hughes [114] calculated the Fermi contact term which yields the hyperfine splitting 
of the ground state of the muonic helium system (a ++ ^“e“) by using a Hylleraas expansion. They 
required hundreds of terms in the expansion because of the slow convergence. Drachman [115] 
noticed that the first term in this expansion gave 99.4% of the contribution, suspecting that a per- 
turbative treatment could be an appropriate way of calculating the Fermi contact term. The large 
ratio of muonic mass to electron mass suggests an adiabatic Born-Oppenheimer approximation. 

The nonrelativistic Hamiltonian of this system is 


H = - 


— V 2 
M x 


-V 


m 



(216) 


x and r are the coordinates of the muon and electron, respectively, relative to the nucleus. The 
reduced masses are M— 201.069 and m= 0.999863 in units of m e . The hyperfine splitting is given 
by 
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Ap = K J <ix(/r'P 2 (x, r)5(x — r), (217) 

where 4/ is the ground state eigenfunction of H , and K— 14196.1 1 MHz. In the Born-oppenheimer 
method, two of the particles are held fixed while we solve for the wave function of the other particle. 
We hold the muon (x) fixed and solve for ^(r) of the electron with x as a parameter 

(-^V 2 - 2 - + V(x, r) - E,)j ^x(r) = 0, (218) 

where U(x, r) = — 2 + | x 2 r | . As indicated in Eq. (106), the solution of Eq. (218) can be written as 


V’x(r) = V>o(r)[l + G(x,r)], 


(219) 


where 


V’o(r) 



(220) 


satisfies the Coulomb problem and G is due to the adiabatic perturbation V(x,r). At this point, 
Drachman expands G in a perturbation series in V and its first order satisfies the equation (108) 


where 


1v 2 Gi -2^- = V- < V >, 
m or 


( 221 ) 


< V >= J dr^ x (r)V^ x (r) = 2[l/a; - m - e~ 2mx (l/x + to)] (222) 

is the expectation value of V . Dalgarno and Lynn [51] have solved an equation similar to Eq. (221) 
but for an electron in the field of two fixed positive charges. Using their solution with suitable 
modification, Drachman obtains 


Gi(x,x) = ^-—x—(2/x)Ei(2x) — (l/x+l)ln(2yx)—e 2 x {i+ 2 x+(l/x+l)[ 21 n( 27 x)-Ei( 2 x)]}, (223) 

^ Li 

where E\ is an exponential integral and Ei is defined in Eq. (114). Now we can solve Eq. (218) for 
E x to obtain E x — —m+ < V > + < VG i> and determine the muonic wave function by solving 
the equation 


= C 224 ) 

Since muon is close to the nucleus due to its large mass, only small values of x are significant and 
the Coulomb term dominates. Therefore ignoring the third term, the solution is hydrogenic: 

$(x) = (8M 3 /7r) l / 2 e ~ 2Mx . (225) 
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Finally, letting $(x, r) = $(x)'0 x (r), Eq. (217) can be written explicitly as 

QO J<T roo 

Av « (Mm) 3 / dxx 2 e~ x / a [l + 2Gi(rax, mx)], (226) 

7T JO 

where a = [4 M + 2m] _1 = 1.24 X 10“ 3 , and the term quadratic in G\ has been dropped for 
consistency, since it is of second order in V. 

Expanding G\ and retaining only the linear term in #, we obtain 


Av (32K /n)(mM) 3 (2a 2 — 12 ma 4 ) 

= (4483.38- 33.36) MHz 

= 4450.02 MHz. (227) 

The quadratic term in G\ contributes 0.689 MHz, while the cubic term is -0.005 MHz. Taking half 
of the quadratic term as an error, the final result is thus Az/=4450.4±0.4 MHz. This agrees with 
the result of Huang and Hughes [114] and is also close to the result obtained by Lakdawala and 
Mohr [117]. 


This is an excellent example of good physics where a rather complicated problem has been 
reduced by Dr. Drachman to a relatively simple problem by realizing that the muon is so close 
to the nucleus that it acts as a hydrogenic system and the rest of the problem is amenable to an 
adiabatic perturbation treatment, where already some available results could be used gainfully to 
obtain fairly accurate results. Understanding the intricacies of a problem and solving it in a simple 
way has been the hallmark of Dr. Drachman ’s research work throughout. 


I wish to thank Keith Feggans for retrieving figures, given in this article, from various published 
papers. 
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